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Abstract. Gamma-ray bursts are the most luminous explosions in the Universe, and 
their origin and mechanism are the focus of intense research and debate. More than 
three decades after their discovery, and after pioneering breakthroughs from space and 
ground experiments, their study is entering a new phase with the recently launched 
Swift satellite. The interplay between these observations and theoretical models of the 
prompt gamma ray burst and its afterglow is reviewed. 
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1. Introduction 

Gamma-ray bursts (GRB) are brief events occurring at an average rate of a few per day 
throughout the universe, which for a brief period of seconds completely flood with their 
radiation an otherwise almost dark gamma-ray sky. While they are on, they outshine 
every other source of gamma-rays in the sky, including the Sun. In fact, they are the most 
concentrated and brightest electromagnetic explosions in the Universe. Until recently, 
they were undetected at any wavelengths other than gamma-rays, which provided poor 
directional information and hence no direct clues about their site of origin. 

This changed in early 1997 when the Beppo-SAX satellite succeeded in detecting 
them in X-rays, which after a delay of some hours yielded sufficiently accurate positions 
for large ground-based telescope follow-up observations. These proved that they were 
at cosmological distances, comparable to those of the most distant galaxies and quasars 
known in the Universe. Since even at these extreme distances (up to Gigaparsecs, or 
~ 10 28 cm) they outshine galaxies and quasars by a very large factor, albeit briefly, their 
energy needs must be far greater. Their electromagnetic energy output during tens of 
seconds is comparable to that of the Sun over ~ few x 10 10 years, the approximate 
age of the universe, or to that of our entire Milky Way over a few years. The current 
interpretation of how this prodigious energy release is produced is that a correspondingly 
large amount of gravitational energy (roughly a solar rest mass) is released in a very short 
time (seconds or less) in a very small region (tens of kilometers or so) by a cataclysmic 
stellar event (the collapse of the core of a massive star, or the subsequent mergers of two 
remnant compact cores). Most of the energy would escape in the first few seconds as 
thermal neutrinos, while another substantial fraction may be emitted as gravitational 
waves. This sudden energy liberation would result in a very high temperature fireball 
expanding at highly relativistic speeds, which undergoes internal dissipation leading 
to gamma-rays, and it would later develop into a blast wave as it decelerates against 
the external medium, producing an afterglow which gets progressively weaker. The 
resulting electromagnetic energy emitted appears to be of the order of a percent or less 
of the total energy output, but even this photon output (in 7-rays) is comparable to the 
total kinetic energy output leading to optical photons by a supernova over weeks. The 
remarkable thing about this theoretical scenario is that it successfully predicts many of 
the observed properties of the bursts. This fireball shock scenario and the blast wave 
model of the ensuing afterglow have been extensively tested against observations, and 
have become the leading paradigms for the current understanding of GRB. 

Historically, GRBs were first discovered in 1967 by the Vela satellites, although they 
were not publicly announced until 1973 [224 . These spacecraft, carrying omnidirectional 
gamma-ray detectors, were flown by the U.S. Department of Defense to monitor for 
nuclear explosions which might violate the Nuclear Test Ban Treaty. When these 
mysterious gamma-ray flashes were first detected, and it was determined that they did 
not come from the Earth's direction, the first suspicion (quickly abandoned) was that 
they might be the product of an advanced extraterrestrial civilization. Soon, however, 
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it was realized that this was a new and extremely puzzling cosmic phenomenon 224J. 
For the next 25 years, only these brief gamma-ray flashes were observed, which could 
be only roughly localized, and which vanished too soon, leaving no traces, or so it 
seemed. Gamma-rays are notoriously hard to focus, so no sharp gamma-ray "images" 
exist to this day: they are just diffuse pin-pricks of gamma-ray light. This mysterious 
phenomenon led to a huge interest and to numerous conferences and publications on 
the subject, as well as to a proliferation of theories. In one famous review article at 
the 1975 Texas Symposium on Relativistic Astrophysics, no fewer than 100 different 
possible theoretical models of GRB were listed |417j . most of which could not be ruled 
out by the observations then available. 

The first significant steps in understanding GRBs started with the 1991 launch 
of the Compton Gamma- Ray Observatory, whose results were summarized in |127j . 
The all-sky survey from the BATSE instrument showed that bursts were isotropically 
distributed, strongly suggesting a cosmological, or possibly an extended galactic 
halo distribution, with essentially zero dipole and quadrupole components |119j . At 
cosmological distances the observed GRB fluxes imply enormous energies, which, from 
the fast time variability, must arise in a small volume in a very short time. This must lead 
to the formation of an e ± — 7 fireball j343t I167| I446j . which will expand relativistically. 
The main difficulty with this scenario was that a smoothly expanding fireball would 
convert most of its energy into kinetic energy of accelerated baryons (rather than 
into photon energy), and would produce a quasi-thermal spectrum, while the typical 
timescales would not explain events much longer than milliseconds. This difficulty was 
addressed by the "fireball shock scenario" 403~| l30lj . based on the realization that shocks 
are likely to arise, e.g. when the fireball ejecta runs into the external medium, after the 
fireball has become optically thin, thus reconverting the expansion kinetic energy into 
non-thermal radiation. The complicated light curves can also be understood, e.g. in 
terms of internal shocks 404, I439[ 122 7j in the outflow itself, before it runs into the 
external medium, caused by velocity variations in the outflow from the source. 

The next major developments came after 1997, when the Italian-Dutch satellite 
Beppo-SAX succeeded in detecting fading X-ray images which, after a delay of 4-6 
hours for processing, led to positions [Hj, allowing follow-ups at optical and other 
wavelengths, e.g. |472j . This paved the way for the measurement of redshift distances, 
the identification of candidate host galaxies, and the confirmation that they were indeed 
at cosmological distances 295| 1100] 1240^ f243j. The detection of other GRB afterglows 
followed in rapid succession, sometimes extending to radio |1321 1134j and over timescales 
of many months |471j . and in a number of cases resulted in the identification of candidate 
host galaxies, e.g. |430| HHJ I341j . etc. The study of afterglows has provided strong 
confirmation for the generic fireball shock model of GRB. This model led to a correct 
prediction [305 , in advance of the observations, of the quantitative nature of afterglows 
at wavelengths longer than 7-rays, which were in substantial agreement with the data 
P%31 147)31 14H51 l47inil5T5] . 

A consolidation of the progress made by Beppo-SAX was made possible through 
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the HETE-2 satellite |195j . after the demise of CGRO and Beppo-SAX. It provided 
a continuing stream of comparable quality afterglow positions, after typical delays of 
hours, and contributed to the characterization of a new class of sources called X-ray 
flashes or XRF |194j resembling softer GRBs, which had been earlier identified with 
Beppo-SAX. It also localized GRB 030329, which resulted in the first unambiguous 
association with a supernova (SN2003dh) 4551 1196j . 

The third wave of significant advances in the field is due to the Swift multi- 
wavelength afterglow satellite, launched in November 2004, which achieved the long- 
awaited goal of accurately localized afterglows starting a minute or so after the burst 
trigger, at gamma-ray, X-ray and optical wavelenghts 46QJ I15fij . This revealed the 
hitherto unexplored afterglow behavior between minutes to hours, enabling a study of 
the transition from the prompt emission and the subsequent long term afterglow, and 
revealing a rich range of X-ray early behavior. It also achieved the long-awaited discovery 
of the afterglows of "short" gamma-ray bursts (whose hard gamma-ray emission is briefer 
than 2 s). It furthermore broke through the symbolic redshift z = 6 barrier, beyond 
which very few objects of any kind have been measured. 

On the theoretical side, a major issue raised by the large redshifts, e.g. |24()[ I244j . is 
that the measured 7-ray fluences (the flux integrated over time) imply a total energy of 
order a solar rest mass, M & c 2 ~ 2 x 10 54 ergs, if it is emitted isotropically. By contrast, 
the total radiant (and the associated kinetic expansion energy) of supernovae (SN), 
which is detected over timescales of weeks to months, is of the order of a thousandth of 
a solar rest mass, 10 51 ergs. A GRB emission which is concentrated in a jet, rather than 
isotropically, alleviates significantly the energy requirements. There is now extensive 
observational evidence for such collimated emission from GRBs, provided by breaks in 
the optical/IR light curves of their afterglows |244t 11401 ^'A ■ The inferred total amount 
of radiant and kinetic energy involved in the explosion is in this case comparable to 
that of supernovae (except that in GRBs the energy is mostly emitted in a jet in 7-rays 
over tens of seconds, whereas in supernovae it is emitted isotropically in the optical over 
weeks). While the luminous (electromagnetic) energy output of a GRB is thus "only" 
of the same order of magnitude as that of supernovae, the explosion is much more 
concentrated, both in time and in direction, so its specific brightness for an observer 
aligned with the jet is many orders of magnitude more intense, and appears at much 
higher characteristic photon energies. Including the collimation correction, the GRB 
electromagnetic emission is energetically quite compatible with an origin in, say, either 
compact mergers of neutron star-neutron star (NS-NS) or black hole-neutron star (BH- 
NS) binaries |343| I105( I331| 1299] . or with a core collapse (hypernova or collapsar) model 
of a massive stellar progenitor |514t I346| 1380} I283| I513j . which would be related to but 
much rarer than core-collapse supernovae. While in both cases the outcome could be, 
at least temporarily, a massive fast-rotating ultra-high magnetic field neutron star (a 
magnetar), the high mass involved is expected to lead inevitably to the formation of a 
central black hole, fed through a brief accretion episode from the surrounding disrupted 
core stellar matter, which provides the energy source for the ejection of relativistic 
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matter responsible for the radiation. 

A stellar origin of GRB leads to two predictions which are similar to those 
for core-collapse supernovae, albeit in so far unobserved aspects. In both GRB 
(whether from compact mergers or from collapsar scenarios) and in core-collapse SN, 
the central material is compressed to nuclear densities and heated to virial temperatures 
characteristically in the multi-MeV range, leading to 5-30 MeV thermal neutrinos. 
And in both cases, the merging or collapsing core material acquires a time-varying 
quadrupole mass moment (which may be smaller in SN not related to GRB), which 
leads to gravitational wave emission. In both GRB and supernovae, the total neutrino 
emission is of the order of a fraction of a solar rest mass, ~ several x 10 53 ergs. The 
gravitational wave emission is of the same order for compact mergers, probably less than 
that for collapsars, and much less in normal core collapse SNe. Experiments currently 
planned or under construction will be able to probe these new channels. 

2. Observational Progress up to 2005 

Before reviewing, in the next section, the latest observational advances achieved with 
Swift, the observational progress made up to that time is briefly surveyed. More 
extensive discussion and references on observations previous to 2004 are, e.g. in 

H23 EZH E23 EZZ1 

The 7-ray phenomenology of GRB was extensively studied and characterized by 
the BATSE instrument on the Compton GRO satellite |127j . The 7-ray spectra are 
non-thermal, typically fitted in the MeV range by broken power-laws whose energy per 
decade peak is in the range 50-500 KeV [TH], sometimes extending to GeV energies 205J. 
GRB appeared to leave no detectable traces at other wavelengths, except in some cases 
briefly in X-rays. The gamma-ray durations range from 10~ 3 s to about 10 3 s, with 
a roughly bimodal distribution of long bursts of t b > 2 s and short bursts of t b < 2s 
[237], and substructure sometimes down to milliseconds. The gamma-ray light curves 
range from smooth, fast-rise and quasi-exponential decay (FREDs), through curves with 
several peaks, to highly variable curves with many peaks |127| I238j (FigureHJ). The pulse 
distribution is complex |369j . and the time histories can provide clues for the geometry 
of the emitting regions |121|. I122j . 

GRB were conclusively shown to be at cosmological distances following Beppo- 
SAX localizations of their X-ray afterglows in 1997 [Hj, followed by optical host galaxy 
identification and redshift determinations |472j . The afterglows decay as a power law in 
time in a manner predicted by pre-existing models [3U5] . softening in time from X-rays 
to optical to radio (e.g. Figure EJ. The energy needed to explain the total (mainly 
gamma-ray) energy fluence can be as large as 10 54 (f2 7 /47r) ergs, where Af2 7 is the solid 
angle into which the gamma-rays are beamed. This is for the highest fluences seen in 
some of the most distant bursts, although for many bursts the energy budget problem is 
not as extreme. If the emission is assumed to be emitted isotropic (isotropic equivalent 
luminosity or energy) this energy ranges up to a solar rest mass in gamma-rays. This 
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Figure 1. Diversity of gamma-ray light curves observed by BATSE 127 



would strain a stellar origin interpretation, since from basic principles and experience 
it is known that, even for the most efficient radiation conversion schemes, a dominant 
fraction of the energy should escape in the form of thermal neutrinos and gravitational 
waves. The energy requirements, however, are much less severe in the case when the 
emission is collimated ( §2.2j) . 

GRB afterglow light curves such as those shown in Figure El have been followed 
up starting several hours after the trigger in X-rays by Beppo-SAX and subsequently 
HETE-2, and in the optical/IR from ground-based telescopes (or in some case with 
HST), and have been explained in terms of forward shock emission (for discussions 
of the pre-Swift data interpretation see, e.g. |471[ I298[ I525[ I377j ). Afterglows have 
been followed up at radio wavelengths in some cases over months, and the analysis and 
interpretation of the radio spectra and light curves |47()^ ESI HHfi] provides important 
clues for the calorimetry and the multi-waveband model fits discussed in $ 15.11 
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Figure 2. The afterglow light curves of GRB 970228 |5T5] . 



2.1. Progenitor candidates 

There is now strong observational evidence [471] that GRB result from a small fraction 
(~ 10~ 6 ) of stars which undergo a catastrophic energy release event toward the end 
of their evolution. For the class of long GRB the candidates are massive stars 
whose core collapses |514[ 1346] I141j to a black hole, either directly or after a brief 
accretion episode, possibly in the course of merging with a companion. This scenario 
is referred to as the collapsar or hypernova scenario, which received strong support 
through the secure spectroscopic detection in some cases of an associated supernova 
event (e.g. |15()| 1455] ll%j : see also §* j2.4| 18. 2j) . For short bursts the most widely 
speculated candidates are mergers of neutron star (NS) binaries or neutron star-black 
hole (BH) binaries (323 1353 QUI HS1 EUl EEEl EUS EHl EMj, which lose orbital 
angular momentum by gravitational wave radiation and undergo a merger. This second 
progenitor scenario has only now begun to be tested thanks to the Swift detection of 
short burst afterglows (see §Bl ©• Both of these progenitor types are expected to 
have as an end result the formation of a few solar mass black hole, surrounded by a 
temporary debris torus whose accretion can provide a sudden release of gravitational 
energy, sufficient to power a burst. An important point is that the overall energetics from 
these various progenitors need not differ by more than about one order of magnitude 
[309] . The duration of the burst in this model is related to the fall-back time of matter to 
form an accretion torus around the BH |141] 1380] or the accretion time of the torus t 332j. 
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Other related scenarios include the formation from a stellar collapse of a fast-rotating 
ultra-high magnetic field neutron star [468 . ffi£\ 13531 I5TJ51 l4T%] . 

2.2. Light curve breaks and jets 

An important subsequent development was the observation, in many of the well-sampled 
afterglows, of a break or steepening of the X-ray and optical light curves [244. 1135| . which 
can be interpreted as being due to the outflow being jet-like and the break occurs when 
the edge of the jet becomes visible, as the jet slows down [408 , 409l l442|l307| . The typical 
(long burst )inferred jet opening angles are 9j ~ 5 — 20 degrees, which reduces the total 
energy requirements from 10 53 — 10 54 erg to ~ 10 51 ergs |135| 1352] . with a dispersion of 
a factor ~ 10. The details are dependent on assumptions about the jet geometry, and 
whether one addresses the gamma-ray energy |135j or the jet kinetic energy [352J, with a 
somewhat larger dispersion in the latter. Variable optical linear polarization is expected 
at the time of a jet break |435[ll65ll415| . which can provide additional constraints, also on 
the jet structure. Light curve break determinations continued through the Beppo-SAX 
and HETE-2 afterglow observation periods, mainly in the optical. Break observations 
in the Swift era are discussed in E J5.51 

2.3. Optical flashes 

Prompt optical flashes (starting within tens of seconds after the gamma-ray trigger) have 
been reported from ground-based small robotic telescopes in a few bursts jHCHl- These 
arise much earlier, are initially brighter and decay more steeply than more ubiquitous 
long-term, slow decaying optical afterglows generally detected since 1997. They are 
also rare: between 1999 and 2004 there were only a handful of prompt optical flashes 
detected with robotic ground telescopes [TOJ . Since the Swift launch, more than twenty 
prompt UVOT (or ground-based robotic) optical flashes have been seen, mostly at times 
starting several hundreds of seconds after the trigger. None have been as bright as the 
first one detected in GRB 990123, except for the notable most distant GRB 050904, 
whose optical brightness is comparable to that of GRB 990123 [54 j. This is discussed 
further in £16.51 

2.4- Association with supernovae 

At least some long GRBs are associated with supernova explosions. The first reported 
example was the GRB980425/SN1998bw association USUI I2H1 IUJ- SN 1998bw was 
a peculiar, energetic Type Ib/c supernova. Using it as a template, other possible 
associations have been claimed through identifying a so-called red supernova bump 
on the optical afterglow light curves of GRB 980326 49 , GRB 970228 ITST]. GRB 
000911(251, GRB 991208[Hnj, GRB 990712[HI|, GRB 011121 50 , GRB 020405 
and GRB 031203 [701 1281] . The first unambiguous supernova signature (SN 2003dh) 
was detected in the z = 0.168 GRB 030329, firmly establishing the GRB-SN 
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associations 053 EH]. Another GRB/SN event, GRB060218/SN2006aj, is discussed 
in and the supernova connection is discussed further in < j8.2l 

2.5. X-ray flashes 

X-ray flashes (XRFs) are a class of bursts whose light curves and spectra resemble 
typical GRB, except for the fact that their spectra are much softer, their spectral peaks 
E p k being typically tens of keV or less |194| [223 . XRFs were first identified with the 
Beppo-SAX satellite, and have been studied in greater detail and numbers with the 
HETE-2 satellite [251 [250 j. Their fluxes and isotropic equivalent luminosities tend to 
be smaller than for GRB, which makes afterglow searches more difficult. Nonetheless, 
several afterglows have been detected, and redshifts have been measured in some of 
them (jSTj, also §$S1EJ- 

2.6. Empirical correlations and distance estimators 

The collimation-corrected total burst energy clustering around 10 51 ergs, while making 
a stellar origin quite plausible, is unfortunately not sufficiently well defined to use 
as standard luminosity candles, whose apparent brightness would provide a distance 
determination. There are other possible distance measures, based on empirical 
correlations between burst observables. One of these is an apparent gamma-ray light 
curve variability correlation with the isotropic equivalent luminosity |12()( I399| . Another 
is the time lag (between higher energy versus lower energy) gamma-rays and the isotropic 
luminosity 336[l2T Hf337j . Attempts at modeling the spectral lags have relied on observer- 
angle dependences of the Doppler boost 321] I432j . In these correlations the isotropic 
equivalent luminosity was used, in the absence of jet signatures, and they must be 
considered tentative for now. However, if confirmed, they could be invaluable for 
independently estimating GRB redshifts. A third one is a correlation between the 
spectral properties and the isotropic luminosity ESI E51 127(ij . These measures can be 
effectively calibrated only for light curves or spectra obtained with the same instrument 
and when redshifts are available. For the above methods the calibration set consists 
typically of a dozen bursts, which is insufficient for being considered reliable. The data 
set is now rapidly increasing with new Swift redshift determinations, but the Swift 
spectra are less well constrained at energies above 150 keV than with Beppo-SAX. 

In recent years, attention has been drawn to a correlation between the photon 
spectral peak energy E p k and the apparent isotropic energy E iso , or the isotropic 
luminosity L iso jSj of the form E p k oc E][^ . For GRB, this correlation has been calibrated 
on a sample consisting of around ten bursts. While strongly suggestive, the sample is 
relatively small and the dispersion is large, so its usefulness as a distance measure is 
precarious. Another suggestive result is that for some X-ray flashes whose redshift has 
been measured the correlation continues to hold [25Qj, even though the peak energy is 
~ 1.5 and the E iso is ~ 3 orders of magnitude lower than for GRB. 
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An interesting development of the E pk — E iso correlation is the proposal of a 
relationship between the collimation-corrected total energy (Ej and the photon spectral 
peak. Modeling the jet break under the assumption that the jet expands into a uniform 
density external medium, this has the form E pk oc E®- 7 for a jet assumed to propagate 
in a homogeneous medium |159| IM| llfilj , and the relation appears to be tighter 
than the previously discussed E pk oc E]^ (Amati) relation. This requires for calibration 
both an observed redshift and a light curve break. If it holds up for a larger calibration 
sample than the current (~ 18 bursts so far), it could be of promise as a cosmological 
tool |25H 1162] . There are however problems to resolve before this becomes competitive 
with SNIa as a cosmological tool |139| . The main ones are the presence of outliers 
[325 , the lack of a large low redshift sample for calibration (since using a high redshift 
sample requires assumptions about the cosmology which it is supposed to test), an 
evaluation of observational biases and selection effects, and the dependence of the results 
on model assumptions about the external medium and jet properties. The latter may 
be circumvented by relying only on observables, e.g. E pea k, fluence (or peak flux) and 
the break time t br [268, I330| . The dispersion, however, remains so far about 1.5-2 times 
larger than for SNIa. 

3. Recent Results from Swift and Follow-up Observations 

Compared to previous missions, the Swift results represent a significant advance on two 
main accounts. First, the sensitivity of the Burst Alert Detector (BAT, in the range 
20-150 keV) is somewhat higher than that of the corresponding instruments in CGRO- 
BATSE, BeppoSAX and HETE-2 [2.T. Second, Swift can slew in less then 100 seconds 
in the direction determined by the BAT instrument, positioning its much higher angular 
resolution X-ray (XRT) and UV-Optical (UVOT) detectors on the burst |15fi] 

As of December 2005, at an average rate of 2 bursts detected per week, over 100 
bursts had been detected by BAT in about a year (compared to 300/year by BATSE; 
note the BAT field of view is 2 sr, versus 4n in BATSE). Of these ~ 100 bursts, 90% 
were detected and followed with the XRT within 350 s from the trigger, and about half 
within 100 s [HE], while ~ 30% were detected also with the UVOT gTT]. Of the total, 
over 23 resulted in redshift determinations. Included in this total sum are nine short 
GRB, of which five had detected X-ray afterglows, three had optical, and one had a 
radio afterglow, with five host galaxy detections and redshift determinations. These 
were the first ever short GRB afterglows detected and followed. 

The new observations bring the total redshift determinations to over 50 since 1997, 
when BeppoSAX enabled the first one. The median redshift of the Swift bursts is z <; 2, 
which is a factor ~ 2 higher than the median of the BeppoSAX and HETE-2 redshifts 
|1U] . This is a statistically significant difference between the Beppo-SAX and Swift 
redshift samples |213^ ITT)] . This may in part be ascribed to the better sensitivity of the 
BAT detector, but mostly to the prompt and accurate positions from XRT and UVOT, 
making possible ground-based detection at a stage when the afterglow is much brighter, 
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by a larger number of robotic and other telescopes. As of January 2006 the highest 
Swift-enabled redshift is that of GRB 050904, obtained with the Subaru telescope, 
z = 6.29 |22()j . and the second highest is GRB 050814 at z — 5.3, whereas the previous 
Beppo-SAX era record was z = 4.5. The relative paucity of UVOT detections versus 
XRT detections may be ascribed in part to this higher median redshift, and in part 
to the higher dust extinction at the implied shorter rest-frame wavelenghts for a given 
observed frequency |411j . although additional effects may be at work too. 

In some of the bursts, both of the "long" (t 7 > 2 s) and "short" (£ 7 < 2 s) categories 
as defined by BATSE, the Swift BAT results show faint soft gamma-ray extensions or 
tails, which extend the duration by a substantial factor beyond what BATSE would 
have detected [156J. A rich trove of information on the burst and afterglow physics has 
come from detailed XRT light curves, starting on average 100 seconds after the trigger, 
together with the corresponding BAT light curves and spectra. This suggests a canonical 
X-ray afterglow picture 338, 1528} IBT] which includes one or more of the following: 

1) an initial steep decay Fx oc t~ ai with a temporal index 3 ^ a± ^ 5, and an 
energy spectrum F u oc v~^ x with energy spectral index 1 < p 1 < 2 (or photon 
number index 2 ^ a + 1 ^ 3), extending up to a time 300 s £ t% ^ 500 s; 

2) a flatter decay portion Fx oc t~ a2 with temporal index 0.2 < a<i ^ 
0.8 and energy index 0.7 ^ /3 2 ~ 1-2, at times 10 3 s ^ t 2 ~ 10 4 s; 

3) a "normal" decay Fx oc t~ a3 with 1.1 < ^ 1.7 and 0.7 ^ (3-s ^ 1.2 (generally 
unchanged the previous stage), up to a time t 3 ~ 10 5 s, or in some cases longer; 

4) In some cases, a steeper decay Fx oc t~ ai with 2 < « 4 < 3, after t 4 ~ 10 5 s; 

5) In about half the afterglows, one or more X-ray flares are observed, sometimes 
starting as early as 100 s after trigger, and sometimes as late as 10 5 s. The energy 
in these flares ranges from a percent up to a value comparable to the prompt emission 
(in GRB 050502b). The rise and decay times of these flares is unusually steep, depending 
on the reference time to, behaving as (t — to) ±afl with 3 ^ otfi ^ 6, and energy indices 
which can be also steeper (e.g. (3fi ^ 1.5) than during the smooth decay portions. 
The flux level after the flare usually decays to the value extrapolated from the value 
before the flare rise (see Figure |3J). The above characteristics are derived mainly from 
long bursts, but interestingly, at least one of the short bursts shows similar features. 
However, the evidence for late time activity is more sketchy in short bursts, so that the 
analogy must be considered with caution. 

An exciting result from Swift was the detection of the long burst GRB 050904, 
which broke through the astrophysically and psychologically significant redshift barrier 
of z ~ 6, which is thought to mark the approximate end of the "dark ages", when re- 
ionization of the intergalactic medium by the first generation of light sources approaches 
completion. This burst was very bright, both in its prompt 7-ray emission (E ljiso ~ 10 54 
erg) and in its X-ray afterglow. Prompt ground-based optical/IR upper limits and 
a J-band detection suggested a photometric redshift z > 6 |188j . and spectroscopic 
confirmation soon followed with the 8.2 m Subaru telescope, giving z = 6.29 220\. 
There are several striking features to this burst. One is the enormous X-ray brightness, 
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Figure 3. Schematic features seen in early X-ray afterglows detected with the Swift 
XRT instrument (e.g. (5281155%] (see text). 



exceeding for a full day the X-ray brightness of the most distant X-ray quasar know to- 
date, SDSS J0130+0524, by up to a factor 10 5 in the first minutes |491| . The implications 
as a tool for probing the intergalactic medium are thought-provoking. Another feature is 
the extremely variable X-ray light curve, showing many large amplitude flares extending 
up to at least a day. A third exciting feature is the report of a brief, very bright IR flash 
[51], comparable in brightness to the famous my ~ 9 optical flash in GRB 990123. 

A third major advance from Swift was the discovery and localization of short GRB 
afterglows. As of December 2005 nine short bursts had been localized by Swift, while in 
the same period HETE-2 discovered two, and one was identified with the IPN network. 
In five of the short bursts, GRB 050509b, 050709, 050724 and 051221a an X-ray afterglow 
was measured and followed up, with GRB 050709, 050724 and 051221a showing also 
an optical afterglow, and 050724 also a radio afterglow, while 040924 had an optical 
afterglow but not an X-ray one |131j . These are the first afterglows detected for short 
bursts. Also, for the first time, host galaxies were identified for these short bursts, which 
in a number of cases are early type (ellipticals) and in other cases are irregular galaxies 
(e.g. |383j . The redshifts of four of them are in the range z ~ 0.15 — 0.5, while another 
one was determined to be z = 0.8 (and less securely, it has been argued that this latter 
may instead be z ~ 1.8 |39j). The median z is ^ 1/3 — 1/2 that of the long bursts. 
There is no evidence for significant star formation in these host environments (except 
for GRB 050709, [130] 1383] . which is compatible with what one expects for an old 
population, such as neutron star mergers or neutron star-black hole mergers, the most 
often discussed progenitor candidates (although it would also be compatible with other 
progenitors involving old compact stars). While the evidence for a neutron star or black 
merger is suggestive, the evidence is not unequivocal. E.g. the observations suggest a 
typical time delay of at least several Gyr between the start formation epoch and the 
explosion of short GRBs 328, 1530j . There are a number of unresolved issues related 
to this (see ^ZJ). The first short burst afterglow followed up by Swift, GRB 050509b, 
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was a rather brief (~ 30 ms), moderate luminosity (L iso ~ 10 50 erg s _1 but low fluence 
(E iso ~ 2 x 10 48 erg) burst with a simple power- law X-ray afterglow which could only 
be followed for several hundreds of seconds |155j . The third one, GRB 050724, was 
brighter, E iso ~ 3 x 10 50 erg, and could be followed in X-rays with Swift for at least 
10 5 s [20], and with Chandra up to 2 x 10 6 s [157] . The remarkable thing about this 
burst's X-ray afterglow is that it shows some similarities to the typical X-ray light curves 
described above for long GRB - except for the lack of a slow-decay phase, and for the 
short prompt emission which places it the the category of short bursts, as well as the 
elliptical host galaxy candidate. It also has also bumps in the X-ray light curve at 100 
s and at 3 x 10 4 s, which resemble some of the long burst X-ray flares and whose origin 
is unclear. The first bump or flare has the same fluence as the prompt emission, while 
the late one has ~ 10% of that. The interpretation of these pose interesting challenges, 
as discussed below and in S 

A fourth exciting result from Swift was the detection of a long burst, GRB 060218, 
which was seen also with the XRT and UVOT instruments [01], and which is associated 
with SN 2006aj [Hill EH01 EMI EH]- The redshift is z = 0.033, and the contribution 
to the optical light curve as well as the spectrum are similar to those of the Ic type 
SN1998bw. The result this is so exciting is that it is the first time that a GRB/SN has 
been observed minutes after the 7-ray trigger at X-ray and UV/optical wavelengths. 
This is discussed further in §8.21 

4. Theoretical Framework 

4-1. The Relativistic Fireball Model 

As discussed in the introduction, the ultimate energy source of GRB is convincingly 
associated with a catastrophic energy release in stellar mass objects. For long bursts, 
this is almost certainly associated with the late stages of the evolution of a massive 
star, namely the collapse of its core |514[ 1346] . which at least in some cases is associated 
with a detectable supernova. For short bursts, it has been long assumed |343t I105j that 
they were associated with compact binary mergers (NS-NS or NS-BH), a view which is 
gaining observational support [4*01 1130j . although the issue cannot be considered settled 
yet. In both cases, the central compact object is likely to be a black hole of several solar 
masses (although it might, temporarily, be a fast rotating high-mass neutron star, which 
eventually must collapse to black hole). In any case, the gravitational energy liberated 
in the collapse or merger involves of order a few solar masses, which is converted into free 
energy on timescales of milliseconds inside a volume of the order of tens of kilometers 
cubed. This prompt energy is then augmented by a comparable amount of energy release 
in a similar or slightly larger volume over a longer timescale of seconds to hundreds of 
seconds, by the continued infall or accretion of gas onto the central object, either from 
the central parts of the massive progenitor star or from the debris of the disrupted 
compact stars which was temporarily held up by its rotation. 
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The principal result of the sudden release of this large gravitational energy (of order 
a solar rest mass) in this compact volume is the conversion of a fraction of that energy 
into neutrinos, initially in thermal equilibrium, and gravitational waves (which are not 
in thermal equilibrium), while a significantly smaller fraction (10~ 2 — 10~ 3 ) goes into a 
high temperature fireball (kT <; MeV) consisting of e ± , 7-rays and baryons. The fireball 
is transparent to the gravitational waves, and beyond several interaction lengths, also 
to the neutrinos. This leads to the prompt emission (on timescales of a few seconds) of 
roughly comparable energy amounts (several xlO 53 ergs) of thermal v e v e with typical 
energies 10-30 MeV, and of gravitational waves mainly near 10 2 — 10 3 Hz. These two, 
by far most dominant, energy forms are so far undetected, and are discussed further in 
A smaller fraction of the liberated energy, or order 10 50 — 10 52 ergs remains trapped 
in a e ± , 7-ray and baryon fireball, which can also contain a comparable (or in some 
scenarios a larger) amount of magnetic field energy. This amount of energy is observed, 
mainly as non-thermal gamma-rays. While smaller than the predicted thermal neutrino 
and gravitational wave fluence, this is nonetheless a formidable electromagnetic energy 
output, much more intense than any other explosive event in the universe. While the 
total energy is comparable to the electromagnetic and kinetic energy of supernovae, the 
difference is that in supernovae the energy is doled out over months, mainly at optical 
wavelengths, while in GRB most of the electromagnetic energy is spilled out in a matter 
of seconds, and mainly at 7-ray wavelengths. 

The leading model for the electromagnetic radiation observed from GRBs is based 
on the relativistic fireball created in the core collapse or merger. The photon luminosity 
inferred from the energies and timescales discussed and from the observations is many 
orders of magnitude larger than the Eddington luminosity Le = AirGMrripC/ ' ot = 
1.25 x 10 38 (M/Mq) erg s -1 , above which radiation pressure exceeds self-gravity, so 
the fireball will expand. The first (thermal) fireball models were assumed to reach 
relativistic expansion velocities [H3J I343[ Ilfi7[ I446j . However, the ultimate expansion 
velocity depends on the baryon load of the fireball |344j . If the fireball energy involved all 
the baryons in the core (solar masses) the expansion would be sub-relativistic. However, 
near the black hole the density is reduced due to accretion and centrifugal forces, it 
is likely that baryons are much depleted in the region where the fireball forms, with a 
tendency to form high entropy (high energy/mass ratio) radiation bubbles. Dynamically 
dominant magnetic fields would also tend to involve fewer baryons. A phenomenological 
argument shows that the expansion must, indeed, be highly relativistic. This is based 
on the fact that most of the GRB spectral energy is observed above 0.5 MeV, so that 
the mean free path for the 77 — > process in an isotropic plasma (an assumption 
appropriate for a sub-relativistically expanding fireball) would be very short. This leads 
to a contradiction, since many bursts show spectra extending above 1 GeV, so the flow 
must be able to avoid degrading these via photon-photon interactions to energies below 
the threshold m e c 2 = 0.511 MeVfTTJT]. To avoid this, it seems inescapable that the flow 
must be expanding with a very high Lorentz factor T, since then the relative angle at 
which the photons collide is less than T _1 and the threshold for the pair production is 
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then diminished. This condition is 

T > 10 2 [(e 7 /10GeV)(ei/MeV)] 1 / 2 , (1) 

in order for photons e 7 to escape annihilation against target photons of energy e t ~ 1 
Mev |297| 1191] . I.e., a relativistically expanding fireball is expected, with bulk Lorentz 
factors T = 100T 2 > 1. 

4-2. Reference frames and timescales in relativistic flows 

The emitting gas is moving relativistically with velocity f3 — v/c — (1 — 1/T 2 ) 1//2 relative 
to a laboratory frame K* which may be taken to be the origin of the explosion or stellar 
frame (which, aside from a cosmological Doppler or redshift factor is the same as the 
Earth frame K of an observer). The lengths, times, thermodynamic and radiation 
quantities of the gas are best evaluated in the gas rest frame (the comoving frame) K', 
and are obtained in the stellar/lab frame through Lorentz transformations. Thus, a 
proper length dr 1 in the comoving frame has a stellar/lab frame length dr* = dr' /V (if 
both ends r*i and r* 2 of the length dr* = r*i — r* 2 in K* are measured at the same 
time so dt* = in K*; i.e. the usual Fitz-Gerald contraction). Similarly, a proper time 
interval dt' in the comoving frame has a duration dt* = dt'T in K* (provided the times 
t*i and t* 2 of dt* = t*i — t* 2 in K* are measured at same positions x*i and x*2 i n K* so 
dx* = 0; the usual time dilation effect). The time needed in the stellar/lab frame K* 
for the gas to move from x*i to x* 2 is the usual dt* = t*i — t* 2 = dr*/ [3c ~ dr*/c. 

When it comes to observations at Earth of radiation emitted from the relativistically 
moving gas, even though the Earth frame K is essentially the same as the K* stellar/lab 
frame, in addition to the above Lorentz transformations one has to consider also the 
classical light travel time delay (Doppler) effect, e.g. |424j . In the observer frame K 
one can use the same spatial coordinates r = r* and dr = dr* as in K*, but the actual 
time of arrival of signals as measured by an observer, which is for brevity denoted just 
t, differs from t* by the above Doppler effect, t ^ t*. Since this observed time t is the 
actual observable, it is customary to describe GRB problems in terms of t (remembering 
it is 7^ t*) and r (which is = r*). Considering a gas which expands radially in a direction 
at an angle cos 8 = fi respect to the observer line of sight, if a first photon is emitted 
when the gas is at the radius r*i = ri (which is at a distance d from the observer) at 
t*i, this photon arrives at the observer at an observer time t\ = t*i — d/c. A second 
photon emitted from a radius r*2 = ^2 at time t*2 will arrive at an observer time 
^2 — t*2 + {d/c — f3fidt*), where dt* = t*i — 1*2- This is illustrated in the source frame in 
Figure EJ For an observer close to the line of sight the observed time difference between 
the arrival of the two photons is 

dt = dt*{l-(3fi) ~ ^(l/2r 2 +fl 2 /2) ~ dr/(2r 2 c)(l+r 2 fl 2 ) ~ dr/(2r 2 c) , (2) 

where we assumed T ^> 1 for an approaching gas (fi = cos 8 > 0) along a radial direction 
well inside the light cone 8 <^ r -1 . 
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Figure 4. Illustration of the emission from spherical relativistic shells in the source 
frame and the relativistic time delay leading to the relation between source frame time 
and observer time. 

While both dt and dt* are in the same reference frame, K = K*, the difference 
is that dt* is the time difference between emission of the two photons, and dt is the 
time difference between reception of the two photons. The general relation between the 
observer frame and comoving frame quantities is given through the Doppler factor £>, 

V=[F(1-^)}- 1 (3) 

where T> ~ 21" for an approaching gas with r ^> 1, /i — > 1 and 9 < (blueshift), or 
T> ~ 1/21" for a receding gas with /i — > — 1 (redshift). Thus, the relation between the 
comoving frame dt' and the observer frame time dt is 

dt = V~ l dt' = T(l - f3n)dt' ~ dt'/2T , (4) 

where an approaching gas is assumed with 9 < (while dt* = dt'T). This is 
illustrated in terms of observer-frame quantities in Figure El Note that in all the 
above transformations we have neglected cosmological effects, which would result in 
multiplying any reception or observer- frame times by an additional factor (1 + z) for 
signals emanating from a source at redshift z. 

The relation between the source frame and observer frame frequency, solid angle, 
specific intensity, temperature, volume, specific emissivity, specific absorption coefficient 
and radial width are obtained in terms of the Doppler factor using relativistic invariants 
HH, v = Vis', dVl = V~ 2 dQ!, I v {y) = V 3 I' u ,(v'), T{y) = VT'(u'), dV = VdV, 
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Figure 5. For a distant observer (located to the right) viewing a shell which expands 
spherically from S with Y = y/l — (v/c) 2 3> 1, the locus of the points from which 
radiation reaches it at a later time t appears as a spheroid (equal arrival time surface) . 
Most of the radiation arrives from the forward (right) hemisphere, which is strongly 
Doppler boosted inside the light cone 1/r (after |4U2| '). The apparent transverse radius 
of the ellipsoid is r± ~ Tci, and its semi-major axis is rii ~ 2Tct, where t is observer 
time. 

jviy) — T^ 2 j' v i{y')i Vviy) — T^~ 1 ^'v'{ v ')i & r = T^~ 1 Sr'. Here \i v = na u (in cm -1 , where n 
is density and a v is absorption cross section), and both v\i v and the optical depth are 
invariants. 

4-3. Relativistic dynamics 

From general considerations, an outflow arising from an initial energy Eq imparted to 
a mass Mo <C Eq/c 2 within a radius tq will lead to an expansion, which due to the 
initial high optical depth can be considered adiabatic. The pressure will be dominated 
by radiation, so the adiabatic index is 7 a = 4/3, and the comoving temperature T' (or 
comoving random Lorentz factor per particle 7') evolves with comoving volume V'as 
T oc V' l ~ la . With a comoving volume V oc r 3 (equation jH]) this means T' oc 7' oc r -1 . 
By conservation of energy, this decrease in internal energy per particle is balanced by 
an increase in its expansion-related energy, i.e. the bulk kinetic energy per particle 
or bulk Lorentz factor T, so that 7r = constant, so that T oc r. This expansion 
occurs at the expense of the comoving frame internal energy. Since the bulk Lorentz 
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factor per particle cannot increase beyond the initial value of random internal energy 
per particle, 70 = rj = E /M c 2 , the bulk Lorentz factor only grows until it reaches 
Tmaz ~ rj = E /M c 2 , which is achieved at a radius r/r ~ 77. Beyond this radius the 
flow begins to coast, with r ~ 77 ~ constant |343| I167| I446| 1345] . 

„/ v f( r / r o), for r/r <r], r < r s ; 

r(r) ~ < (5) 
[ 77 , for r/r > 77, r >r s 

which defines a saturation radius r s ~ rorj beyond which the Lorentz factor has 
saturated. Another way to understand the initial acceleration |3UUj is that initially, 
at r = ro, the gas particles have a bulk Lorentz factor T ~ 1 and have an isotropical 
distribution of velocities with random Lorentz factors 7 ~ 77 = E /M c 2 . As the 
particles expand outward, when they have reached a radius r their velocity vectors 
will confined inside an angle (r/Yo) -1 of the radial direction. A transformation to a 
comoving frame moving radially with a bulk Lorentz factor T(r) ~ r/r is needed for 
the velocity distribution to be isotropic in the comoving frame, as it should be. 

As particles initially contained inside r move outwards with velocity vectors which 
are increasingly radial, they form a radially expanding shell whose lab-frame width 
is initially 5r ~ 5r Q ~ r . The radial velocity spread is (c — v)/c = 1 — /3 ~ T~ 2 , 
which causes a gradual spread of the lab-frame radial width 5r/r ~ 5v/v ~ T~ 2 . For 
typical values of ro ~ 10 6 — 10 7 cm and 77 < 10 3 this is negligible until well beyond the 
saturation radius, and a noticeable departure from the approximately constant width 
5r ~ r Q starts to become appreciable only for radii in excess of a spreading radius rs 
where 5r ~ rbvjc ~ rsi]~ 2 > r . The laboratory frame width is therefore |3UUj 

r r, /F 2l / 6r ° ' fOT r ~ r ^ 

6r ~ max [5r , r/T \ ~ J (6) 

I r/1 , for r > r 5 . 

where the spreading radius rs ~ Sr^r] 2 is a factor 77 larger than the saturation radius 
r s ~ 5r ?7. 

The comoving radial width 5r' is related to the lab width 5r through 8r' = 5rT> ~ 
5rr. Hence 

{5r T ~ r for r < r s ; 

5r r] for r s < r < r 5 ; (7) 

r/77 for r > rs , 

Since the dimensions transverse to the motion are invariant, the comoving volume 
is V oc r 2 5r', which behaves as 300 

{Anr 3 for r < r s ; 

Air7]5r r 2 for r s < r < rs; (8) 

47T77~ 1 r 3 for r > r^; 

and the comoving particle density n' oc V'~ l . For an adiabatic expansion (valid for the 
high initial optical depths) and a relativistic gas polytropic index 4/3 (valid as long as 
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the pressure is dominated by radiation), one has 

i/Q [{& r o/r) for r < r s ; 

E' T { V \ ' 

— = — = ( J ^ < (5r D /r s )(r s /r) 2/3 = r/" 1/3 (5r /r) 2/3 for r s < r < r 5 ; (9) 

-fO / „,„ 

[ (5r /r s )(r s /r s ) 1 (r 5 /r) = rj 1 (5r Q /r) for r > r 5 ; 

where E', T', p', V are comoving internal energy, temperature, density and volume. 

The above equations refer to the release of an energy Eq and mass Mq corresponding 
to rj — Eo/Mqc 2 , originating inside a region of dimension 5tq ~ tq. This mass and 
energy leaves that original region in a lab-frame (or observer frame) light-crossing 
time 5to ~ Sr /c. For typical core collapse or compact merger stellar scenarios, the 
energy release volume is of the order of several Schwarzschild radii of the ensuing 
black hole (BH), few times 2GM BH /c 2 with M BH > 2M Q , say r ~ 10 7 mi cm, where 
mi = M BH /10M Q , with a light-crossing timescale t ~ r /c ~ 3 x 10~ 4 m! s. This is 
of the order of the dynamical (Kepler) timescale near the last stable circular orbit in a 
temporary accretion disk feeding the newly formed black hole (or near the light cylinder 
of an initial fast-rotating magnetar or neutron star, before it collapses to a black hole). 

4-4- Optical Depth and Photosphere 

As shown by 343 , 1167[ I446j a relativistically expanding fireball initially has e ± pairs 
in equilibrium which dominate the scattering optical depth, but the pairs fall out of 
equilibrium and recombine below a comoving temperature T" ~ 17 keV, and thereafter 
only a residual freeze-out density of pairs remains, which for r\ not too large (in practice 
rj ^ 10 5 i.e. baryon loads not too small) is much less than the density of "baryonic" 
electrons associated with the protons, n e = n p . For a typical burst conditions the initial 
black-body temperature Tq at ro ~ 10 7 cm is a few MeV, and pair recombination occurs 
at radii below the saturation radius. The scattering optical depth of a minishell (and 
of the whole outflow) is still large at this radius, due to the baryonic electrons. For a 
minishell of initial width 5ro the optical depth varies as [300J 

r T = (M a T /47rr 2 m p ) = T (5r /r) 2 (10) 

where r = (E ctt / ' ^r\m v c 2 r}) , for 5r Q ^ r/T 2 or r ^ r$- Assuming a burst with 
total energy E = 10 52 E 52 and total duration t gr b divided into minishells of duration 
8t = 3 x 10~ 4 s, each of energy 10 47 " 5 Kt7.5 erg, these becomes optically thin at 

r t = 7- 1/2 r = 3 x lOVlW*) 1 ^ cm ; (11) 

where henceforth the notation Q x (where a; is a number) indicates the quantity Q in 
units of 10 x times its c.g.s. units. 

For bursts of some substantial duration, e.g. an outflow duration t gr b = 10 s as 
above, at any instant different parts of the flow have different densities, and are above or 
below the saturation radius, so a continuous outflow picture is more appropriate 344J. 
In this "wind" regime one defines the dimensionless entropy as rj = L/Mc 2 , and instead 
of integral conservation laws one uses the relativistic fluid differential equations. The 
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Lorentz factor again grows linearly and saturates at the same radius r s = r 7] (equation 
jHJ), where r = 5r Q is the minimum variability radius, and the adiabatic behavior of 
equation (jHJ) is the same for the temperature, etc. The particle density follows from the 
mass conservation equation, 

n ' p = (M/4irr 2 m p cT) = (L/A7vr 2 m p c 3 7]T) (12) 

and the optical depth is r T (r) = / r °° n' e a T [(l - /3)/(l + f3)} 1/2 dr' ~ n' e a T (r/2T) which 
yields for the global photosphere 344J 

r ph ~ (Ma T /%irm p cT 2 ) = (La T /8iTm p c 3 rir 2 ) ~ 6 x 10 n L 5 i% 3 cm. (13) 

The comoving temperature of the flow behaves as T' oc r _1 , r~ 2 / 3 and the observer- 
frame temperature T = T'Y is T ~ T , T ~ T (r/r s )~ 2 / 3 for r < r s , r > r s (equations 

ED- 

The radiation escaping from a radius r (e.g. the photospheric radius r p h) which is 
released at the same stellar frame time t* would arrive at the observer only from within 
angles inside the light cone, 9 < The observer-frame time delay between light 

coming from central line of sight and the edges of the light cone (the so-called angular 
time |167j ) is 

t ang ~ t ~ (r/c)(l - (3) ~ r/2cT 2 . (14) 

This is because the 'edge' of the light cone corresponds to an angle of 1/r from the 
line of sight, and therefore ct ang ~ r(l — cos^) ~ r(l — (3) ~ r/2T 2 , since at 6* = 1/r, 
cos^ ~ (3. This time is the same as the observer- frame time of equation (J2J). Note that if 
the outflow duration t gr b is shorter than t ang of equation (flljl. the latter is the observed 
duration of the photospheric radiation (due to the angular time delay). Otherwise, for 
tgrb > t ang , the photospheric radiation is expected for a lab- frame duration t grb . 

4-5. Thermal vs. Dissipative Fireballs and Shocks 

The spectrum of the photosphere would be expected to be a black-body |343[I167[ I446J. 
at most modified by comptonization at the higher energy part of the spectrum. However, 
the observed 7-ray spectrum observed is generally a broken power law, i.e., highly non- 
thermal. In addition, a greater problem is that the expansion would lead to a conversion 
of internal energy into kinetic energy of expansion, so even after the fireball becomes 
optically thin, it would be highly inefficient, most of the energy being in the kinetic 
energy of the associated protons, rather than in photons. For a photosphere occurring 
at r < r s , which requires high values of 7], the radiative luminosity in the observer frame 
is undiminished, since E' rad oc r -1 but T oc r so E rad ~ constant, or Lph oc r 2 T 2 T' 4 oc 
constant, since T' oc r _1 . However for the more moderate values of rj the photosphere 
occurs above the saturation radius, and whereas the kinetic energy of the baryons is 
constant E kin ~ E ~ constant the radiation energy drops as E rad oc (r/r s )~ 2 / 3 , or 
L ph ~ L (r ph /r s )- 2 / 3 [Ml HUH! ■ 

A natural way to achieve a non-thermal spectrum in an energetically efficient 
manner is by having the kinetic energy of the flow re-converted into random energy 
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via shocks, after the flow has become optically thin 0D31 EQU EEB C2I3 EHI 033 • Such 
shocks will be collisionless, i.e. mediated by chaotic electric and magnetic fields rather 
than by binary particle interactions, as known from interplanetary experiments and as 
inferred in supernova remnants and in active galactic nuclei (AGN) jets . As in these 
well studied sources, these shocks can be expected to accelerate particles via the Fermi 
process to ultra- relativistic energies [I3J 12 H08[ I265| I45U] , and the relativistic electron 
component can produce non-thermal radiation via the synchrotron and inverse Compton 
(IC) processes. A shock is essentially unavoidable as the fireball runs into the external 
medium, producing a blast wave. The external medium may be the interstellar medium 
(ISM), or the pre-ejected stellar wind from the progenitor before the collapse. For an 
outflow of total energy E and terminal coasting bulk Lorentz factor T = rj expanding 
in an external medium of average particle density n ex t, the external shock becomes 
important at a deceleration radius r^ec for which E = (47r / 3)r dec n m p c 2 ri 2 [403], 

r dec ~ {ZEvl^n^m^ri 1 ) 1 ' 3 ~ 5.5 x W^E^n; 1 / 3 ^ 3 cm . (15) 

At this radius the initial bulk Lorentz factor has decreased to approximately half its 
original value, as the fireball ejecta is decelerated by the swept-up external matter. 
The amount of external matter swept at this time is a fraction rj^ 1 of the ejecta 
mass Mo, M ext ~ Mo/77 (in contrast to the sub- relativistic supernova expansion, where 
deceleration occurs when this fraction is ~ 1). 

The light travel time difference between a photon originating from r = and a 
photon originating from matter which has moved to a radius r with a Lorentz factor 
T is At ~ (r/c)(l — /3) ~ r/2cr 2 |4U2j . and the emission from a photosphere or from 
a shock emission region at radius r moving at constant T is also received from within 
the causal light cone angle on an observer angular timescale t ~ r/2cr 2 |403t 1301]. 
For an explosion which is impulsive (i.e. essentially instantaneous as far as observed 
relativistic time delays) a similarity solution of the relativistic flow equations shows that 
the bulk of the ejected matter at a radius r is mainly concentrated inside a region of 
width Ar ~ r/2r 2 |44[ I45j . The time delay between radiation along the central line of 
sight originating from the back and front edges of this shell also arrive with a similar 
time delay t ~ r/2cr. Thus, the timescale over which the deceleration is observed to 
occur is generally 

t dec ~ W(2cr 2 ) ~ lO^K) 1 / 3 ^ 3 s, (16) 

and this is the observer timescale over which the external shock radiation is detected. 
This is provided that the explosion can be taken to be impulsive, which can be defined as 
the outflow having a source-frame (and observer frame) duration t gr b < td ec (see however 
Variability on timescales shorter than tdec may occur on the cooling timescale or 
on the dynamic timescale for inhomogeneities in the external medium, but is not ideal 
for reproducing highly variable profiles [439 j . and may therefore be applicable to the class 
of long, smooth bursts. However, it can reproduce bursts with several peaks |3~53] . and 
if the external medium is extremely lumpy (An /n 10 5 — 10 6 ) it might also describe 
spiky GRB light curves jHSj. 
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Before the ejecta runs into the external medium, "internal shocks" can also occur as 
faster portions of the ejecta overtake slower ones, leading to pp collisions and 7r° decay 
gamma-rays |387| and to fast time- varying MeV gamma-rays |4()4j . The latter can be 
interpreted as the main burst itself. If the outflow is described by an energy outflow 
rate L Q and a mass loss rate M Q = dM /dt starting at a lower radius rj, maintained over 
a time T, then the dimensionless entropy is r] = L /M c 2 , and the behavior is similar 
to that in the impulsive case, T oc r and comoving temperature T" oc r^ 1 , followed by 
saturation T max ~ 77 at the radius r/r D ~ 77 |344| . For variations of the output energy 
or mass loss of order unity, the ejected shells of different Lorentz factors A77 ~ r] are 
initially separated by ct v (where t v < T are the typical variations in the energy at ri), 
and they catch up with each other at an internal shock (or dissipation) radius 

rdis ~ ct v rf ~ 3 x 10 14 ^ i0 ?72 cm > ( 17 ) 

The time variability should reflect the variability of the central engine, which might 
be expected e.g. from accretion disk intermittency, flares, etc. |331j . The radiation 
from the disk or flares, however, cannot be observed directly, since it occurs well below 
the scattering photosphere of the outflow and the variability of the photons below it is 
washed out |4U4j . The comoving Thomson optical depth is tt = n' e aTr/T, and above 
the saturation radius r s = r r\ where T = r], the radius of the photosphere (tt = 1), is 
given from equation (|13|) as 

!1.2 x lO 12 ^!^ 3 cm, for r > r s ; 

1.2 x lO^Lsfrof^ 1/3 cm for r <r s . ' 

The location of this baryonic photosphere defines a critical dimensionless entropy 
i]* = 562(L 51 /r 7) 1//4 above (below) which the photosphere occurs below (above) 
the saturation radius |31()j . In order for internal shocks to occur above the wind 
photosphere and above the saturation radius (so that most of the energy does not 
come out in the photospheric quasi-thermal radiation component) one needs to have 
3.3 x lO^L^roj/t^o) 1 / 5 < r] < 5.62 x 10 2 (L 5 iAo, 7 ) 1/4 . The radial variation of the bulk 
Lorentz factor and the location of the various characteristic radii discussed above is 
shown in Figure El 

Such internal shock models have the advantage |4U4j that they allow an arbitrarily 
complicated light curve, the shortest variation timescale t V)m i n ^ 10~ 4 s being limited 
only by the dynamic timescale at tq ~ ct V)Tn i n ~ 10 7 ro,7 cm, where the energy input 
may be expected to vary chaotically, while the total duration is tg^ ^> t v . Such 
internal shocks have been shown explicitly to reproduce (and be required by) some of 
the more complicated light curves [4"3U| 1227} l358j (see however |96| B25J). The gamma- 
ray emission of GRB from internal shocks radiating via a synchrotron and/or inverse 
Compton mechanism reproduces the general features of the gamma-ray observations 
[1381 145 fj . There remain, however, questions concerning the low energy (20 — 50 keV) 
spectral slopes for some bursts (see §sec:spectrum). Alternatively, the main 7-ray bursts 
could be (at least in part) due to the early part of the external shock [403, 96J. Issues 
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arise with the radiation efficiency, which for internal shocks, is estimated to be moderate 
in the bolometric sense (5-20%), higher values 30 — 50%) being obtained if the shells 
have widely differing Lorentz factors |451[ I3*T| E32 , although in this case one might 
expect large variations in the spectral peak energy E peak between spikes in the same 
burst, which is problematic. The total efficiency is substantially affected by inverse 
Compton losses |362[ 1374] . The efficiency for emitting in the BATSE range is typically 
low ~ 1 — 5%, both when the MeV break is due to synchrotron |245H45l|ll85| and when 
it is due to inverse Compton 354J. 




Figure 6. Jet Lorentz factor schematic behavior and examples of nominal locations 
of the saturation radius r s , photospheric radius r p ^, internal shock (or magnetic 
dissipation) radius Ti S and external shock r es . The photosphere produces thermal 
7-rays, the internal shock/dissipation region produces the non-thermal 7-rays, the 
external shock region produces the afterglow. 



4-6. Duration, reverse shocks, thin and thick shells 

In the following discussion we assume for simplicity a uniform external medium. For 
a baryonic outflow such as we have been considering, the timescale to ~ r o/ c ~ ms 
represents a minimum variability timescale in the energy-mass outflow. (Note, however, 
if the gamma-ray emission arises from local dissipation events, such as e.g. magnetic 
reconnection in a Poynting flux dominated outflow, the minimum timescales could be 
smaller than the timescales of the central source variations). On the other hand, the 
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total duration t grb of the outflow, during which the central engine keeps pouring out 
energy and matter, is likely to be substantially longer than the minimum variability 
timescale to- The temporary accretion disk must have an outer radius larger than ro, 
and a total accretion (or jet energization) time t gr b ^> to ( or the magnetar has a spin- 
down time tg r b to)- Thus, in general the total lab-frame width of the outflow ejecta 
will be A « ctgrb, which may be viewed as composed of many radial minishells whose 
individual widths are 5r ~ cSto or larger. While the saturation radius is still r s ~ 5r ?7 
where 5tq ~ ro corresponds to the shortest variability time (and the smallest minishells 
coast after this r s ), the entirety of the ejecta reaches coasting speed only after its leading 
edge has moved to a larger radius ?v ~ Ar?, and the ejecta as whole starts to spread 
at a larger radius ~ Arj 2 (even though individual minishells of initial width 5r <C A 
start to spread individually at the smaller radius rs ~ Sr r) 2 ). 

In general, whatever the duration t gr b of the outflow, one expects the external shock 
to have both a forward shock (blast wave) component propagating into the external 
medium, and a reverse shock propagating back into the ejecta |3()lj . The forward 
shock and the reverse shock start forming as soon as the outflow starts, although their 
radiation is initially weak and increases progressively. The forward shock is highly 
relativistic, r ~ rj from the very beginning, but the reverse shock starts initially as a 
sub-relativistic sound wave (relative to the contact discontinuity or shock frame) and 
becomes progressively stronger as more external matter is swept up. (This describes 
the more frequently encountered "thin shell" case, see below; the reverse shock becomes 
stronger with time only if the external density profile is shallower than r~ 2 , whereas the 
reverse shock strength is constant for an r~ 2 profile at r < r^). 

For an impulsive regime outflow, where t gr b < td ec , i-e. when the outflow time is 
shorter than the time-delayed duration of the external shock when it starts to decelerate, 
equation \Ug\, this deceleration time can be taken to be the observable duration of the 
peak emission from the external shock. Thereafter the expansion goes into a self-similar 
expansion with T oc r -3 / 2 [44, 403j. In this case, td ec is also the observer time at which 
the reverse shock finishes crossing the ejecta, and at that time the reverse shock Lorentz 
factor f r relative to the contact discontinuity frame has become marginally relativistic, 
r r ~ 1, while relative to the external gas or the observer, the reverse shocked gas is still 
moving at almost the same speed as the forward shocked gas |3()HI3()4| . One consequence 
of this is that while the forward shocked protons have highly relativistic random Lorentz 
factors, those in the reverse shock are marginally relativistic, and consequently the 
electrons in the forward shock are much more relativistic than those in the reverse 
shock, leading to a much softer (optical) spectrum of the reverse shocks |302[ 1305] (see 

However, when the outflow time t gr b exceeds the deceleration time t dec of equation 
(llfij) . the eternal shock dynamics is different |438j . In this case there is an initial 
intermediate regime T oc r -1 / 2 (obtained, for a constant external density p and a 
constant kinetic luminosity L at t < t gr b from momentum balance in the shock frame, 
Lj (r 2 r 2 cx pr 2 ), and the transition to a self-similar expansion T oc r~ 3 / 2 j44[ I433j occurs 
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at the observer time t gr b, instead of at tdec- Thus, the observer time for the transition 
to the self-similar expansion is 

T = max[V, 10(E 53 /n ) 1/3 rill /3 s] (19) 
This defines a critical initial Lorentz factor To ~ of the burst by setting T = t gr b in 
place of tdec in equation (JTfij) . 

T BM ~ 300(SB3/no) 1/8 (T/10 s)" 3 / 8 . (20) 

For i] < Tbm, T = tdec we have the usual "thin shell" case, where deceleration and 
transition to the self-similar expansion occurs at the usual r^ ec , td ec , and at this time the 
reverse shock has crossed the ejecta and is marginally relativistic. For 77 > Tbm we have 
a "thick shell" case, where deceleration and transition to the self-similar regime occurs 
at T = t gr b and Tbm ~ 2> C TT BM , when T ~ Tbm- In this 77 > T^m case, the reverse 
shock becomes relativistic, and by the time it has crossed the ejecta (at time T = t gr b) 
the reverse shock Lorentz factor in the contact discontinuity frame is T r ~ i]/2Tbm ^ 1> 
and the forward shock Lorentz factor at this time is T ~ Tbm- 

4-7. Spectrum of the Prompt GRB Emission 

The prompt emission observed from classical GRB (as opposed to XRFs or SGRs) has 
most of its energy concentrated in the gamma-ray energy range 0.1-2 MeV. The generic 
phenomenological photon spectrum is a broken power law with a break energy in the 
above range, and power law extensions down into the X-ray, and up into the 100 MeV 
to GeV ranges (although a substantial fraction of GRB have soft X-ray excesses above 
this, and some are classified as X-ray rich (XRR) [381], a classification intermediate 
between XRF and GRB). For classical GRB the photon energy flux F E oc E-P has 
typical indices below and above the typical observed break energy Eb r ~ 0.2 MeV 
of fl\ ~ and /3 2 ~ 1 JH]- (Pre-BATSE analyses sometimes approximated this as a 
bremsstrahlung-like spectrum with an exponential cutoff at Eb r , but BATSE showed that 
generally the extension above the break is a power law). A synchrotron interpretation 
is thus natural, as has been argued e.g. since the earliest external shock synchrotron 
models were formulated. 

The simplest synchrotron shock model starts from the conditions behind the 
relativistic forward shock or blast wave |403[ 1301 . The post-shock particle and internal 
energy density follow from the relativistic strong shock transition relations [44J, 

n 2 = (4r 2 i + 3)ni ^ 4r 2 i7ii, 

62 = (r 2 i - l)n 2 m p c 2 ~ T 2 in 2 m p c 2 ~ AT^nxnipC 2 , (21) 

where it is assumed that the upstream material is cold. Here n is number density 
and e is internal energy density, both measured in the comoving frames of the fluids, 
r 2 i — r is the relative Lorentz factor between the fluids 2 (shocked, downstream) and 1 
(unshocked, upstream), and the Lorentz factor of the shock front itself is = \/2T 2 i, 
valid for T 2 i 3> 1. For internal shocks the jump conditions can be taken approximately 
the same, but replacing T 21 by a lower relative Lorentz factor r r ~ 1. 
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The typical proton crossing a strong shock front with a relative bulk Lorentz factor 
r 2 i acquires (in the comoving frame) an internal energy characterized by a random 
(comoving) Lorentz factor 7 p m ~ V |3()1| . The comoving magnetic field behind the 
shock can build up due to turbulent dynamo effects behind the shocks |301|. |302j (as 
also inferred in supernova remnant shocks). More recently, the Weibel instability has 
been studied in this context 1291] I335[ 1294} [152 . While the efficiency of this process 
remains under debate, one can parametrize the resulting magnetic field as having an 
energy density behind the shock which is a fraction e# of the equipartition value relative 
to the proton random energy density behind the shock, B' ~ \52ii6Bn ex ('y' p — l)m p c 2 ] 1/,2 r, 
where the post-shock proton comoving internal energy is (7' — l)m p c 2 ~ 1 (or ~ T) 
for internal (external) shocks |301| H04J. Scattering of electrons (and protons) by 
magnetic irregularities upstream and downstream can lead to a Fermi acceleration 
process resulting in a relativistic power law distribution of energies ^(7) oc 7~ p with 
p > 2. It should be stressed that although the essential features of this process are 
thought to be largely correct, and it is widely used for explaining supernova remnant, 
AGN and other non-thermal source radiation spectra, the details are only sketchily 
understood, [HI EJ QUB1 EH3 EZD USD]. (Possible difficulties with the simplest version 
of Fermi acceleration and alternative possibilities were discussed, e.g. in |2*3*1 I192[ ITlOj). 
The starting minimum (comoving) Lorentz factor of the thermal electrons injected into 
the acceleration process, 7e )m would in principle be the same as for the protons, T, 
(they experience the same velocity difference), hence both before and after acceleration 
they would have ~ (m e /m p ) less energy than the protons. However, the shocks being 
collisionless, i.e. mediated by chaotic electric and magnetic fields, can redistribute 
the proton energy between the electrons and protons, up to some fraction e e of the 
thermal energy equipartition value with the protons, so 7 6im ~ e e {m p /m e )Y \602\ l304j. 
If only a fraction Q e < 1 of all the shocked thermal electrons is able to achieve this e e 
initial equipartition value to be injected into the acceleration process, then the initial 
minimum electron random comoving Lorentz factor is 7 m ~ (e e /( e )(m p /m e )T jHO], where 
henceforth we ignore the subscript e in 7e im . More accurately, integrating over the power 
law distribution, one has 7 m = g(p)(m p /m e )(e e /( e )T ~ 310[gQo)/(l/6)](e e /C e )r, where 
g(j?) = (p — 2)/(p — 1). The observer frame synchrotron spectral peak is 

v m ~ r(3/8ir)(eB'/m e c)j 2 m ~ 2 x 10 6 B'j 2 m T Hz , (22) 

and the optically thin synchrotron spectrum is 



v 



1/3 for v < u, 



F^{'_ ln _^ (23) 



v 



-(p-l)/2 for v > 



m 



assuming that the radiative losses are small (adiabatic regime). For the forward external 
shock at deceleration, typical values are, e.g. B' ~ S0(eB,~i n exY^ 2r ]2.5 G, T ~ 77 ~ 3xl0 2 , 
7 m ~ 10 5 (e e _i/C e ,-i)77 2 .5 and v m ~ 2 x 10 20 (e ei _i/C e ,-i) 2 (e i j i _in e:c ) 1 / 2 ^ 5 Hz, while for 
internal shocks typical values are, e.g. B' ~ 3 x 10 5 (eB i _in' 13 ) 1 / 2 r re / i o G, T re i ~ ir re ; j0 , 
7m~3x lO^e^i/Ce.-Or^o and i/ m ~2x 10 19 (e ei _ 1 /^eta ei _ 1 ) 2 (e Bi _in'13) 1/2 r 3 r7 2 .5 
Hz. For the prompt emission, the high energy slope 02 = (p — l)/2 is close to the mean 
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high energy slope of the Band fit, while the low energy slope can easily approach /3x ~ 
considering observations from, e.g., a range of B' values (a similar explanation as for 
the flattening of the low energy synchrotron slope in flat spectrum radio-quasars). The 
basic synchrotron spectrum is modified at low energies by synchrotron self-absorption 
[302] I304| 12 17\ 1171] , where it makes the spectrum steeper (F u ~ v 2 for an absorption 
frequency v a < u m ). It is also modified at high energies, due to inverse Compton effects 
[M EMI EH EH EH EH!, extending into the GeV range. 

The synchrotron interpretation of the GRB radiation is the most straightforward. 
However, a number of effects can modify the simple synchrotron spectrum. One is 
that the cooling could be rapid, i.e. when the comoving synchrotron cooling time 
t' sy = 9mgC 5 /4e 4 -B /2 7 e ) ~ 7 x 10 8 / B' 2 ^ e s is less than the comoving dynamic time 
t' dyn ~ r/2cT, the electrons cool down to 7 C = 67im e c/aTB' 2 t' dyn and the spectrum above 
u c ~ r(3/87r)(eS'/m e c)7 2 is F u oc v~ 1/2 |4^H|[TH4] . Also, the distribution of observed low 
energy spectral indices (3\ (where F v oc u 131 below the spectral peak) has a mean value 
(3\ ~ 0, but for a fraction of bursts this slope reaches positive values (3i > 1/3 which 
are incompatible with a simple synchrotron interpretation |381j . Possible explanations 
include synchrotron self-absorption in the X-ray |172j or in the optical range up-scattered 
to X-rays [354J, low-pitch angle scattering or jitter radiation |292| 1293] . observational 
selection biases |274j and/or time-dependent acceleration and radiation [275J, where 
low-pitch angle diffusion can also explain high energy indices steeper than predicted by 
isotropic scattering. Other models invoke a photospheric component and pair formation 
[310 , see below. 

There has been extensive work indicating that the apparent clustering of the break 
energy of prompt GRB spectra in the 50-500 keV range may be real |381j . rather than 
due to observational selection effects |375j . I.e. the question is, if this is a real clustering, 
what is the physical reason for it. (Note, however, that if X-ray flashes, or XRF, 
discussed below, form a continuum with GRB, then this clustering stretches out to 
much lower energies; at the moment, however, the number of XRFs with known break 
energies is small). Since the synchrotron peak frequency observed is directly dependent 
on the bulk Lorentz factor, which may be random, the question arises whether this peak 
is indeed due to synchrotron, or to some other effect. An alternative is to attribute a 
preferred peak to a black-body at the comoving pair recombination temperature in 
the fireball photosphere |106j . In this case a steep low energy spectral slope is due 
to the Rayleigh- Jeans part of the photosphere, and the high energy power law spectra 
and GeV emission require a separate explanation. For such photospheres to occur at 
the pair recombination temperature in the accelerating regime requires an extremely 
low baryon load. For very large baryon loads, a related explanation has been invoked 
[465 , considering scattering of photospheric photons off MHD turbulence in the coasting 
portion of the outflow, which up-scatters the adiabatically cooled photons up to the 
observed break energy. 

Pair formation can become important [404, 362, 374J in internal shocks or 
dissipation regions occurring at small radii, since a high comoving luminosity implies a 
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large comoving compactness parameter 

£' = n\a T r dis /T ~ {aLa T /4nr dis m e c 3 T 3 ) > 1 , (24) 

where a ^ 1 is the luminosity fraction above the electron rest mass. Pair-breakdown 
may cause a continuous rather then an abrupt heating and lead to a self-regulating 
moderate optical thickness pair plasma at sub-relativistic temperature, suggesting a 
comptonized spectrum |164| . Copious pair formation in internal shocks may in fact 
extend the photosphere beyond the baryonic photosphere value ()18|) . A generic model 
has been proposed |310[ I312( I390| I426~( 142 7\ 13*94] which includes the emission of a thermal 
photosphere as well as a non-thermal component from internal shocks outside of it, 
subject to pair breakdown, which can produce both steep low energy spectra, preferred 
breaks and a power law at high energies. A moderate to high scattering depth can lead 
to a Compton equilibrium which gives spectral peaks in the right energy range |364[ f365j. 
An important aspect is that Compton equilibrium of internal shock electrons or pairs 
with photospheric photons lead to a high radiative efficiency, as well as to spectra with 
a break at the right preferred energy and steep low energy slopes |406( 1366( 1367 . It also 
leads to possible physical explanations for the Amati jH] or Ghirlanda |16Uj relations 
between spectral peak energy and burst fluence |406| [464 . 

4-8. Alternative Prompt Emission Models 

There are several alternative models for the prompt GRB emission, which so far have 
not found wide use for explaining the observations. The most plausible of these, despite 
the technical difficulties which impair its applicability, considers the main 7-ray burst 
emission to arise from magnetic reconnection or dissipation processes, if the ejecta is 
highly magnetized or Poynting dominated ~EM EES EOS EOEl CHI E2H1 MM- The 
central engine could also in principle be a temporary highly magnetized neutron star 
or magnetar |5U8j . These scenarios would lead to alternative dissipation radii, instead 
of equation (fT7]) . where reconnection leads to particle acceleration, and a high radiative 
efficiency is in principle conceivable due to the very high magnetic field. An external 
shock would follow after this, whose radius in the "thin shell" limit would be again given 
by equation ("I"""") , with a standard forward blast wave but no (or a weaker) reverse shock 
304[ I305j . due to the very high Alfven speed in the ejecta. For a long duration outflow, 
however, the dynamics and the deceleration radius would be similar to the "thick shell" 
case of ^4.61 i.e. the case with a relativistic reverse shock [279J. Following the claim 
of an observed high gamma-ray polarization in the burst GRB 021206 [22!, there was 
increased attention on such models for some time (e.g. |279p . and on whether the usual 
baryonic (i.e. sub-dominant magnetization) jets might also be able to produce such high 
polarization "HU H221 [123 E221 12HD1 12S3 [Ml EH- The issue may remain unresolved, 
as the observational analysis appears to be inconclusive 423, 12U 1511] . 

Other alternative models include different central energy sources such as strange 
stars ( [6*o"ll"""""."l01|l348| ) an d charged black hole electric discharges |421j . while retaining 
essentially similar fireball shock scenarios. A model unifying SGR, XRF and GRB 
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|115| Hlb] postulates a very thin (10~ 4 rad) precessing, long-lived magnetized jet. This 
requires a separate explanation for the light-curve ("jet") breaks, and the interaction 
during precession with the massive stellar progenitor is unclear. Another speculative 
radiation scenario considers non-fluid ejecta in the form of discrete "bullets" [193J, or 
"cannon-balls" ejected at relativistic velocities, which assume no collective interactions 
(i.e. no collisionless shocks) and instead rely on particle-particle interactions, and 
produce prompt emission by blue-shifted bremsstrahlung and produce afterglows by IC 
scattering progenitor or ambient photons [ZZHHH] • The predictions are similar to those of 
the standard fluid jet with shocks or dissipation. However, the basic ansatz of coherent 
bullet formation, acceleration to relativistic velocities and their survival against plasma 
instabilities is an unanswered issue in this model. It is also farther from astrophysical 
experience, whereas other well-observed systems such as AGN jets, which are known to 
be fluid (as is almost everything else in astrophysics at high energy per particle values) 
involve dynamical and radiation physics concepts which are quite plausibly extended to 
the GRB context. Fluid or plasma GRB outflow and jet models are better supported by 
theoretical work and simulations, and are so far not only compatible with observations 
but have produced predictions borne out by observations. Nonetheless, even in this 
standard scenario, the models remain largely phenomenological. The detailed nature of 
the underlying central engine and progenitor are poorly known, and the micro-physics 
of particle acceleration, magnetic field amplification in shocks and/or reconnection or 
dissipation is not well understood, and the radiation mechanisms are, at least for the 
prompt emission, subject of discussion. 

5. Afterglow Radiation Models 

5.1. The standard model 

The external shock starts to develop as soon as the ejecta expands into the external 
medium. As the ejecta plows ahead, it sweeps up an increasing amount of external 
matter, and the bolometric luminosity of the shock increases as L oc t 2 (equating in 
the contact discontinuity frame the kinetic flux L/^nr 2 T 2 to the external ram pressure 
p ex tT 2 while r ~ T = rj ~ constant, r oc 2T 2 ct oc t |403j). The luminosity peaks after 
T has dropped to about half its initial value, at a radius r^ ec at an observer time td ec 
given by equations (jl5flfijh Thereafter, as more matter is swept up, the bulk Lorentz 
factor and the radius vary as [403, 1356j as 

r oc r" 3/2 oc r 3/8 , roct 1/4 (adiabatic), 



or in general T oc r~ 9 oc t~9/( 1+2 s) ; r oc t 1 ^ 14 " 2 ^ with g = (3,3/2) for the radiative 
(adiabatic) regime. In the adiabatic case the radiative cooling time, e.g. synchrotron, 
is longer than the observer-frame dynamical time t ~ r/2cT 2 , so the energy is 
approximately conserved E = (47r /3)r 3 nQm p c 2 T 2 ~ constant (c.f. equation ^3]), while 
in the radiative case the cooling time is shorter than the dynamic time and momentum 



r oc r" 3 oc r 3/7 



r oc t 1/r 




(25) 
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is conserved (as in the snow-plow phase of supernova remnants), n r 3 T ~ constant. 
Thus, after the external shock luminosity peaks, one expects the bolometric luminosity 
to decay as L oc in the adiabatic regime |403j or steeper in the radiative regime, in a 
gradual fading. The observed time-radius relation is more generally t ~ r/KcT 2 , where 
K = 2 in the constant T regime, and K = 4 in the self-similar (BM) regime |497| 1433] . 

The spectrum of radiation is likely to be due to synchrotron radiation, whose peak 
frequency in the observer frame is v m oc ^B'j^, and both the comoving field B' and 
electron Lorentz factor 7 e are likely to be proportional to 7 |3()lj . This implies that as 7 
decreases, so will v mi and the radiation will move to longer wavelengths. Consequences 
of this are the expectation that the burst would leave a radio remnant |347j after some 
weeks, and before that an optical [218J transient. The observation of linear polarization 
at the few percent level observed in a number of optical or IR afterglows (e.g. |471j ) 
supports the paradigm of synchrotron emission as the dominant emission mechanism in 
the afterglow. 

The first self-consistent afterglow calculations [305 took into account both the 
dynamical evolution and its interplay with the relativistic particle acceleration and a 
specific relativistically beamed radiation mechanism resulted in quantitative predictions 
for the entire spectral evolution, going through the X-ray, optical and radio range. For 
a spherical fireball advancing into an approximately smooth external environment, the 
bulk Lorentz factor decreases as in inverse power of the time (asymptotically t~ 3 / 8 in the 
adiabatic limit), and the accelerated electron minimum random Lorentz factor and the 
turbulent magnetic field also decrease as inverse power laws in time. The synchrotron 
peak energy corresponding to the time- dependent minimum Lorentz factor and magnetic 
field then moves to softer energies as r 3 / 2 . These can be generalized in a straightforward 
manner when in the radiative regime, or in presence of density gradients, etc.. The radio 
spectrum is initially expected to be self-absorbed, and becomes optically thin after ~ 
hours. For times beyond ~ 10 minutes, the dominant radiation is from the forward 
shock, for which the flux at a given frequency and the synchrotron peak frequency 
decay as [305J 




as long as the expansion is relativistic. This is referred to as the "standard" (adiabatic) 
model, where g = 3/2 in T oc r~ 9 and (3 = d log F u /d log v is the photon spectral energy 
flux slope. More generally |3()7j the relativistic forward shock flux and frequency peak 
are given by 



where g = (3/2,3) for the adiabatic (radiative) regime. The transition to the non- 
relativistic expansion regime has been discussed, e.g. by |515|IHT|l273j . A reverse shock 
component is also expected |302| 13051 144 1[ [307J, with high initial optical brightness 
but much faster decay rate than the forward shock, see §5.2j) . Remarkably, the simple 
"standard" model where reverse shock effects are ignored is a good approximation for 
modeling observations starting a few hours after the trigger, as during 1997-1998. 




F v oc 




(26) 



F v OC t [3-2 fl (l-2/3)]/(l+2fl) and Vm K r 4 9 /(l+2 s ) 



(27) 



CONTENTS 



32 



The afterglow spectrum at a given instant of time depends on the flux observed at 
different frequencies from electrons with (comoving) energy 7 e m c c 2 and bulk Lorentz 
factor T, whose observed peak frequency is v — T^ 2 (eB' /2irm e c). Three critical 
frequencies are defined by the three characteristic electron energies. These are u m (the 
"peak" or injection frequency corresponding to 7 m ), v c (the cooling frequency), and vm 
(the maximum synchrotron frequency). There is one more frequency, u a , corresponding 
to the synchrotron self-absorption at lower frequencies. For a given behavior of T with 
r or t (e.g. adiabatic, T oc r~ 3 / 2 ) and values of the isotropic equivalent kinetic energy of 
the explosion, of the electron index (e.g. p = 2.2) and the efficiency factors e e , ( e , e^, 
one can obtain the time dependence of the characteristic observer-frame frequencies, 
including also a cosmological redshift factor z 



v m = (6 x 10 15 Hz) (1 + z) l l 2 g{p) 2 {e e /Q 2 efEli 2 t- d i/2 (28) 

v c = (9 x 10 12 Hz) (1 + z)- l/2 t B V2 n- l E^ l2 r d 1/2 (29) 

v a = (2 x 10 9 Hz) (1 + z)-\e e /Q- l eT n^E 1 ^ (30) 

F Ujmax = (20 mJy) (1 + z)e^V/ 2 £ 52 ^ 2 28 , (31) 

where tj = (t/day) and g(p) = (p — 2)/(p — 1). The final GRB afterglow synchrotron 
spectrum is a four-segment broken power law 440, 30HJ I175| 1525] separated by the 
typical frequencies v a , v m , and u c (Figure |7j). Depending on the order between v m and 
u c , there are two types of spectra 440J. For v m < z/ c , called the "slow cooling case", the 
spectrum is 

K/^m) 1/3 (zV^a) 2 V < V a 



A vjmax 



F 

A u,max 



[ylv m )-^l 2 v m <v<v c (32) 

{v c lv m Y^I\vlv c Y^ 2 V C <V< V M 

For v m > v c , called the "fast cooling case", the spectrum is 

KM) V: VK) 2 » < "a 

(v/v c y/ 3 v a <v<v c 
(v/vc)- 1 / 2 v c <v<v m 
{y m lv c )~ xl2 {ylv m )~ vl2 v m < v < v M 

A useful tabulation of the temporal indices a and spectral indices j3 is given in 
Table 1 of |525j , corresponding to the various forward shock spectral regimes of equations 
(J3*2|) . (}3*3*j) . for a homogeneous or a wind external medium. In the above, the normalization 
Fornax is obtained by multiplying the total number of radiating electrons 4nr 3 ni/3 by the 
peak flux from a single electron |440j . which is only a function of B and is independent 
of the energy (7 e ) of the electron |440^ I51()j . There are more complicated regimes for 
various cases of self-absorption |172j . e.g. there can also be an intermediate fast cooling 
optically thick power law segment of the synchrotron spectrum where F v oc v 11 ^ 8 . 

The predictions of the fireball shock afterglow model j305] were made in advance 
of the first X-ray detections by Beppo-SAX [71] allowing subsequent follow-ups 
|472| 1295^ H33j over different wavelengths, which showed a good agreement with the 
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Figure 7. Fast cooling and slow cooling synchrotron spectra |44()| 

standard model, e.g. 4831 15151 1463} I495|. I496|. |400j. The comparison of increasingly 
sophisticated versions of this theoretical model (e.g. [4401 151U1 137(j| W7\ 1175] ) against 
an increasingly detailed array of observations (e.g. as summarized in |471j ) has provided 
confirmation of this generic fireball shock model of GRB afterglows. 

A snapshot spectrum of the standard model at any given time consists generally of 
three or four segment power law with two or three breaks, such as those shown in Figure 
13 (More rarely, a five segment power law spectrum may also be expected |172j ). The 
observations (e.g. |471j ) are compatible with an electron spectral index p ~ 2.2 — 2.5, 
which is typical of shock acceleration, e.g. |495[I44()( 15*10] . etc. As the remnant expands 
the photon spectrum moves to lower frequencies, and the flux in a given band decays 
as a power law in time, whose index can change as the characteristic frequencies move 
through it. Snapshot spectra have been deduced by extrapolating measurements at 
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different wavelengths and times, and assuming spherical symmetry and using the model 
time dependences [496, 1510j . fits were obtained for the different physical parameters of 
the burst and environment, e.g. the total energy E, the magnetic and electron-proton 
coupling parameters e# and e e and the external density n . These lead to typical values 
n Q ~ KT 2 - 10 cm" 3 , e B ~ 10~ 2 , e e ~ 0.1 - 0.5 and E ~ 10 52 - 10 54 ergs (if spherical; 
but see ^5.5|) . 

5.2. Prompt Flashes and Reverse Shocks 

An interesting development was the observation [3] of a prompt and extremely bright 
{m v ~ 9) optical flash in the burst GRB 990123, the first data point for which was 
at 15 seconds after the GRB started (while the gamma-rays were still going on). This 
observation was followed by a small number of other prompt optical flashes, generally not 
as bright. A prompt multi- wavelength flash, contemporaneous with the 7-ray emission 
and reaching such optical magnitude levels is an expected consequence of the reverse 
component of external shocks |302t I305| 144 1[ I307j . Generally the reverse shock can 
expected to be mildly relativistic (thin shell case; see, however, below). In this case the 
thermal Lorentz factor of the reverse electrons is roughly Y e ~ e e m p /m e (whereas in the 
forward shock, the thermal Lorentz factor of the electrons is % ~ e e Tm p /m e . In this 
case the reverse electrons radiate much softer radiation than the forward shock electrons. 
This follows also from the fact that the reverse shock has a similar total energy as the 
forward shock, but consists of V times more electrons, hence the energy per electron is 
1/r times smaller 305J. In general, since the pressure (and hence the magnetic energy 
density) is the same in the forward and reverse shocked regions, one has the following 
relations between forward and reverse shock radiation properties |437j : 1) The peak flux 
of the reverse shock, at any time, is larger by a factor of V than that of the forward shock, 
Flmax = ^Flmax'i 2) The typical frequency of the minimal electron in the reverse shock 
is smaller by a factor of T 2 , v T m = z/^/r 2 ; 3) The cooling frequency of the reverse and 
forward shock are equal, v r c = v{ = v c (under the assumption that es is the same in the 
forward and reverse shocked gas; this might not be true if the ejecta carries a significant 
magnetic field from the source); 4) Generally (also in refreshed shocks) v T ^ < v%f and 
v a < v c- The self-absorption frequency of the reverse shock is larger than that of the 
forward shock. The characteristic frequencies and flux temporal slopes for a standard 
afterglow are given by the case (r) with s = in Table 1 below. 

The prompt optical flashes, starting with GRB 990123, have been generally 
interpreted |44H I307[ I327j as the radiation from a reverse (external) shock, although a 
prompt optical flash could be expected from either an internal shock or the reverse part 
of the external shock, or both |305j . The decay rate of the optical flux from reverse 
shocks is much faster (and that of internal shocks is faster still) than that of forward 
shocks, so the emission of the latter dominate after tens of minutes |169j . Such bright 
prompt flashes, however, appear to be relatively rare. Other early optical flashes, e.g. 
in GRB 021004, GRB 021211, GRB 041219a, GRB 050904 are also consistent with the 
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Table 1. Temporal exponents of the peak frequency v m , the maximum flux F^ ml 
the cooling frequency v c and the flux in a given bandwidth F v , for the forward 
(f) and reverse (r) shocks, calculated both in the adiabatic regime v m < v < v c 
(F v oc F Vm (y m j v)P cx t~ a v~P ', where /3 = (p — l)/2), and in the cooling regime 
v c < v m < v (F v oc (vc/vmY^ivm/v) 13 cx t~ a v~^ where (3 = p/2). For s = 1 this gives 
the usual (i.e. without "refreshment") forward and reverse shock behavior [4371 14(j£] . 



reverse shock interpretation (UHl E2B EE2E1 1123 HUE HH E03 - After the launch of 
Swift, new prompt optical observations with robotic telescopes have greatly added to 
the phenomenology of prompt flashes (see ^21 

5.3. Dependence on external density, injection variability and anisotropy 

If the external medium is inhomogeneous, e.g. n ex t oc r~~ k , the energy conservation 
condition is r 2 r 3 ~ fc ~ constant, so T oc t 1 ^ 4- ^, r oc t-( 3 - k )/( 8 ~ 2k ) ( which changes the 
temporal decay rates 308 . This might occur if the external medium is a stellar wind 
from the evolved progenitor star of a long burst, e.g. n ext oc r~ 2 , such light curves fitting 
some bursts better with this hypothesis [HSII267J. 

Another departure from a simple injection approximation is one where E (or L ) 
and To are not a simple a delta function or top hat functions. An example is if the 
mass and energy injected during the burst duration t gr b (say tens of seconds) obeys 
M(> 7) oc 7~ s , E(> 7) oc 7 1 ~ s , i.e. more energy emitted with lower Lorentz factors at 
later times, but still shorter than the gamma-ray pulse duration |4()5[ 1437] . The ejecta 
dynamics becomes 

r( r ) oc r -(3-*0/(l+s) fx £-(3-fc)/(7+s-2fc) ? r ^ £ (l+s)/(7+s-2fc) _ 

This can drastically change the temporal decay rate, extending the afterglow lifetime 
in the relativistic regime. If can provide an explanation for shallower decay rates, 
if the progressively slower ejecta arrives continuously, re-energizing the external 
shocks ("refreshed" shocks) on timescales comparable to the afterglow time scale 
[405] I247| IH2 E37J. While observational motivations for this were present already in 
the Beppo-SAX €?rct ; ctS discussed in the above references, this mechanism has been 
invoked more recently in order to explain the Swift prompt X-ray afterglow shallow 
decays (see ffii.2j) . When the distribution of Y is discontinuous, it can also explain a 
sudden increase in the flux, leading to bumps in the light curve. After the onset, the 
non-standard decay rates for the forward and reverse shock are tabulated for different 
cases |437j in Table [T] 

Other types of non-standard decay can occur if the outflow has a transverse 9 
dependent gradient in the energy or Lorentz factor, e.g. as some power law E oc 
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9~ a , r oc 9~ b |308j . Expressions for the temporal decay index a(/3, s, d, a, 6, ..) in F u oc t a 
are given by |308| &37J, which now depend also on s, d, a, b, etc. (and not just on j3 
as in the standard relation of equ. pS]) . The result is that the decay can be flatter (or 
steeper, depending on s, d, etc)) than the simple standard a = (3/2)/?. Such non- 
uniform outflows have been considered more recently in the context of jet breaks based 
on structured jets ( fl5.5|) . 

Evidence for departures from the simple standard model was present even before 
the new Swift observations, by e.g. sharp rises or humps in the light curves followed 
by a renewed decay, as in GRB 970508 [363 , 1378] . or shallower than usual light curve 
decays. Time-dependent model fits [359 to the X-ray, optical and radio light curves 
of GRB 970228 and GRB 970508 indicated that in order to explain the humps, a non- 
uniform injection or an anisotropic outflow is required. Another example is the well- 
studied wiggly optical light curve of GRB 030329, for which refreshed shocks provide 
the likeliest explanation |176j . Other ways to get light curve bumps which are not too 
steep after ~ hours to days is with micro-lensing |154[ I173j . late injection |522| I21()j . 
or inverse Compton effects [4361 15231 I190j . The changes in the shock physics and 
dynamics in highly magnetized or Poynting dominated outflows were discussed, e.g. 
in |468| I465| 1306} 1174} II 76} I279| 1527] . More examples and references to departures 
from the standard model are discussed, e.g. in j471[ 1525] . Departures from spherical 
symmetry and jet effects are discussed in the next two subsections. 

5.4- Equal arrival time surface and limb brightening effect 

As illustrated in FigureEl for a distant observer the photons from a spherically expanding 
shell are received from an equal-arrival time surface which is an ellipsoid (if V = 
constant). The photons arriving from the line of sight originated at larger radii than 
photons arriving from the light-cone at 9 ~ T. At smaller radii the outflow had a higher 
magnetic field and higher density, so the radiation from the 1/T edge is harder and more 
intense. Thus an interesting effect, which arises even in spherical outflows, is that the 
effective emitting region seen by the observer resembles a ring |496| 1168] 1357} 14*3*4*] H70j . 
This limb brightening effect is different in the different power law segments of the 
spectrum. When one considers the change in T due to deceleration, the ellipsoid is 
changed into an egg shape, which is similarly limb-brightened. This effect is thought 
to be implicated in giving rise to the radio diffractive scintillation pattern seen in 
several afterglows, since this requires the emitting source to be of small dimensions 
(the ring width), e.g. in GRB 970508 |498j . This provided an important observational 
check, giving a direct confirmation of the relativistic source expansion and a direct 
determination of the (expected) source size |498|l219j . The above treatments were based 
on the simple asymptotic scaling behavior for the Lorentz factor T ~ constant at r < r^ec 
and T oc r -3 / 2 (r oc r~ 3 ) at r > r^ ec for the adiabatic (fully radiative) cases (i J4.5|) . More 
exact treatments are possible [4*T] 14*2*] based on following analytically and numerically 
the detailed dynamical evolution equations for the Lorentz factor through and beyond 
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the transition between pre-deceleration and post-deceleration. The shape of the equi- 
temporal surfaces is modified, and the expected light curves will be correspondingly 
changed. The exact afterglow behavior will depend on the unknown external medium 
density and on whether and what kind of continued of continued energy injection into 
the shock occurs, which introduces an additional layer of parameters to be fitted. 

5.5. Jets 

The spherical assumption is valid even when considering a relativistic outflow collimated 
within some jet of solid angle fij < Air, provided the observer line of sight is inside 

— 1/2 

this angle, and T > |3()()j . so the light-cone is inside the jet boundary (causally 

disconnected) and the observer is unaware of what is outside the jet. However, as the 
ejecta is decelerated, the Lorentz factor eventually drops below this value, and a change 
is expected in the dynamics and the light curves |4U8| $i09\. It is thought that this is 
what gives rise to the achromatic optical light curve breaks seen in many afterglows 

The jet opening angle can be obtained form the observer time tj at which the flux 
F v decay rate achromatically changes to a steeper value, assuming that this corresponds 
to the causal (light-cone) angle T(t) -1 having become comparable to (and later larger 
than) the jet half-angle 0j [408J. Assuming a standard adiabatic dynamics and a uniform 
external medium, the jet opening half-angle is 

9 j ~ 5deg tfjE^n^^/O^iil + z}/2)-^ (35) 

where -E53 is the isotropic equivalent gamma-ray energy in ergs, tjj = tj /day and 
r/ 7 is radiative efficiency |135j . The degree of steepening of the observed flux light 
curve can be estimated by considering that while the causal angle is smaller than 
the jet opening angle, the effective transverse area from which radiation is received 
is A ~ r\ ~ (r/r) 2 oc t 2 T 2 , whereas after the causal angle becomes larger than the jet 
angle, the area is A ~ r 2 6 2 . Thus the flux after the break, for an adiabatic behavior 
T oc t~ 3 / 8 (valid if there is no sideways expansion) is steeper by a factor cx T 2 oc t -3 / 4 
307], a value in broad agreement with observed breaks. After this time, if the jet 
collimation is simply ballistic (i.e. not due to magnetic or other dynamical effects) the 
jet can start expanding sideways at the comoving (relativistic) speed of sound, leading 
to a different decay T oc t~ 1/2 and F v oc t~ p oc t~ 2 [4713] . 

A collimated outflow greatly alleviates the energy requirements of GRB. If the 
burst energy were emitted isotropically, the energy required spreads over many orders 
of magnitude, E ltiso ~ 10 51 — 10 54 erg |243j . However, taking into account the jet 
interpretation of light curve breaks in optical afterglows |352| 11351 13f>3| the spread in 
the total 7-ray energy is reduced to one order of magnitude, around a less demanding 
mean value of E ljtot ~ 1.3 x 10 51 erg [51j . This is not significantly larger than the kinetic 
energies in core-collapse supernovae, but the photons are concentrated in the gamma-ray 
range, and the outflow is substantially more collimated than in the SN case. Radiative 
inefficiencies and the additional energy which must be associated with the proton and 
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magnetic field components increase this value (e.g. the ?7 7 factor in equation [55]). but 
this energy is still well within the theoretical energetics ^ io 53 - 5 — 10 54 erg achievable in 
either NS-NS, NS-BH mergers [HOE! or in collapsar models [SHI EH1 EHH! using MHD 
extraction of the spin energy of a disrupted torus and/or a central fast spinning BH. It 
is worth noting that jets do not invalidate the usefulness of spherical snapshot spectral 
fits, since the latter constrain only the energy per solid angle 309J . 

Equation ()35|) assumes a uniform external medium, which fits most afterglows, but 
in some cases a wind-like external medium (n ext oc r~ 2 ) is preferred |351[ EEl P67J. 
For an external medium varying as n ext = Ar~ k one can show that the the Lorentz 
factor initially evolves as V oc (E/v4) 1 / 2 r~( 3_fc)//2 oc ( J E/A) 1 /(8-2fe) t -(3-fe)/(8-2fc) ; and the 
causality (or jet break) condition Y ~ 9j leads to a relation between the observed light 
curve break time tj and the inferred collimation angle 9j which is different from equation 
(J35J), namely^- oc (E/A)- 1 /(8-2fe)(t j /[i + 2 ])(3-fc)/(8-2/=) K (£/4)-i/4( f j /[l + z ]) 1 /4 > where 

the last part is for k = 2. Another argument indicating that the medium in the vicinity 
of at least some long-GRB afterglows is not stratified, e.g. as r -2 , is the observation of 
a sharp jet-break in the optical afterglow lightcurves (as, e.g. in GRB 990510, 000301c, 
990123). As pointed out by [248 , relativistic jets propagating in a wind-like external 
medium are expected to give rise to a very gradual and shallow break in the afterglow 
light curve. 

The discussion above also makes the simplifying assumption of a uniform jet 
(uniform energy and Lorentz factor inside the jet opening angle, or top-hat jet model). 
In this case the correlation between the inverse beaming factor / 6 ~ = (9 2 /2)~ 1 (or 
observationally, the jet break time from which 9j is derived) and the isotropic equivalent 
energy or fluence E ljiso is interpreted as due to a distribution of jet angles, larger angles 
leading to lower E Jiiso , according to E lyiso oc 9~ 2 . There is, however, an equally plausible 
interpretation for this correlation, namely that one could have a universal jet profile 
such that the energy per unit solid angle dE 1 /dVL oc 9~ 2 , where 9 is the angle measured 
from the axis of symmetry |414t I524j . (To avoid a singularity, one can assume this 
law to be valid outside some small core solid angle). This model also explains the 
|1351 1352j correlation, the different E iso would be due to the observer being at different 
angles relative to the jet axis. This hypothesis has been tested in a variety of ways 
1861 1324( I249| IT7H] . Attempts to extend the universal 9~ 2 jet structure to include X-ray 
flashes ( ^2.5|) . together with use of the Amati relation between the spectral peak energy 
Ep ea k and E 1 ^ so ( §2.fij) . leads to the conclusion that a uniform top-hat model is preferred 
over a universal 9~ 2 jet model 250 . Uniform jets seen off-axis have also been considered 
as models for XRF in a unified scheme, e.g. |517[ I181j . On the other hand, another 
type of universal jet profile with a Gaussian shape |520| IS3] appears to satisfy both the 
jet break- E lyiso and E pea k — E 7t i so correlations for both GRB and XRFs. More extensive 
discussion of this is in [525J. 

The uniform and structured jets are expected to produce achromatic breaks in 
the light curves, at least for wavebands not too widely separated. However, in some 
bursts there have been indications of different light curve break times for widely 
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separated wavebands, e.g. GRB 030329, suggesting different beam opening angles for 
the optical/X-ray and the radio components 37 . Such two-component jets could arise 
naturally in the collapsar model, e.g. with a narrow, high Lorentz factor central jet 
producing 7, X-ray and optical radiation, and a wider slower outflow, e.g. involving 
more baryon-rich portions of the envelope producing radio radiation [392 . A wider 
component may also be connected to a neutron- rich part of the outflow |370j . More 
recent discussions of possible chromatic breaks are in |113t I36lj . 

6. Current Theoretical Issues in Afterglow Models 

The afterglow is generally assumed to become important after the time when the 
self-similar Y oc r~ 3 / 2 behavior starts. From equation (fT9"l for the deceleration time 
tdec ~ ijdecl^cY 2 ) and taking into account the gradual transition to the self-similar 
regime [235]. this is approximately 

t ag ~ T ~ Max[V(l + z) , 15(E 53 /n ) 1 /3r ? - 5 8 / 3 [(l + z )j2\ s] , (36) 

where t gr b is the duration of the outflow, i.e. an upper limit for the duration of the 
prompt 7-ray emission, and a cosmological time dilation factor is included. (Note 
that in some bursts the 7-rays could continue in the self-similar phase). The afterglow 
emission from the forward and the reverse shock emission starts immediately after the 
ejecta starts to sweeps up matter, but it does not peak (and become dominant over the 
prompt emission or its decaying tail) until the time ~ t ag , marking the beginning of the 
self-similar blast wave regime. 

Denoting the frequency and time dependence of the afterglow spectral energy flux 
as F v (t) oc u~^t~ a , the late X-ray afterglow phases (3) and (4) of §H]seen by Swift are 
similar to those known previously from Beppo-SAX (the theoretical understanding of 
which is discussed in and in [525]). The "normal" decay phase (3), with temporal 
decay indices a ~ 1.1 — 1.5 and spectral energy indices (3 ~ 0.7 — 1.0, is what is generally 
expected from the evolution of the forward shock in the Blandford-McKee self-similar 
late time regime, under the assumption of synchrotron emission. 

6.1. Early steep decay 

Among the new afterglow features detected by Swift (see Figure EJ), the steep initial 
decay phase F u oc t~ 3 — i~ 5 in X-rays of the long GRB afterglows is one of the most 
striking. There are several possible mechanisms which could cause this. The most 
obvious first guess would be to attribute it to the cooling following cessation of the 
prompt emission (internal shocks or dissipation). If the comoving magnetic field in the 
emission region is random [or transverse], the flux per unit frequency along the line of 
sight in a given energy band, as a function of the electron energy index p, decays as 
F v oc t~ a with a = —2p [(1 — 3p)/2] in the slow cooling regime, where (3 = (p — l)/2, 
and it decays as a = —2(1 +p), [—(2 — 3p)/2] in the fast cooling regime where (3 = p/2, 
i.e. for the standard p = 2.5 this would be a = —5, [—3.25] in the slow cooling or 
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a = —7, [—2.75] in the fast cooling regime, for random [transverse] fields |307| . In some 
bursts this may be the explanation, but in others the time and spectral indices do not 
correspond well. In addition, if the flux along the line of sight decays as steeply as 
above, the observed flux would be dominated by the so-called high latitude emission, 
which is discussed next. 

At present, the most widely considered explanation for the fast decay, both of the 
initial phase (1) and of the steep flares, attributes it to the off-axis emission from regions 
at 9 > r _1 (the curvature effect, or high latitude emission [226j . In this case, after the 
line of sight gamma-rays have ceased, the off-axis emission observed from 9 > is 
(T#) -6 smaller than that from the line of sight. Integrating over the equal arrival time 
region, this flux ratio becomes oc (T#)~ 4 . Since the emission from 9 arrives (T#) 2 later 
than from 9 = 0, the observer sees the flux falling as F v oc t~ 2 , if the flux were frequency 
independent. For a source-frame flux oc z/ -/3 , the observed flux per unit frequency varies 
then as 

F u <x(t-t y 2 -? (37) 

i.e. a = 2 + (3. This "high latitude" radiation, which for observers outside the line 
cone at 9 > would appear as prompt 7-ray emission from dissipation at radius r, 
appears to observers along the line of sight (inside the light cone) to arrive delayed by 
t ~ (r9 2 /2c) relative to the trigger time, and its spectrum is softened by the Doppler 
factor oc into the X-ray observer band. For the initial prompt decay, the onset of 
the afterglow (e.g. phases 2 or 3), which also come from the line of sight, may overlap 
in time with the delayed high latitude emission. In equation (J37]l to can be taken as the 
trigger time, or some value comparable or less than equation (jHEJ)- This can be used to 
constrain the prompt emission radius |257j . When td ec < T, the emission can have an 
admixture of high latitude and afterglow, and this can lead to decay rates intermediate 
between the two |34Uj . Values of t closer to the onset of the decay also lead to steeper 
slopes. It is possible to identify for various bursts values of to near the rising part of 
the last spike in the prompt emission which satisfy the subsequent steep decay slope 
270 . Structured jets, when viewed on-beam produce essentially the same slopes as 
homogeneous jets, while off-beam observing can lead to shallower slopes |l()3j . For the 
flares, if their origin is assumed to be internal (e.g. some form of late internal shock or 
dissipation) the value of to is just before the flare, e.g the observer time at which the 
internal dissipation starts to be observable [526 . This interpretation appears, so far, 
compatible with most of the Swift afterglows |528[ I338| |360 . 

Alternatively, the initial fast decay may be due to the emission of a cocoon of 
exhaust gas |368j . where the temporal and spectral index are explained through an 
approximately power-law behavior of escape times and spectral modification of multiply 
scattered photons. The fast decay may also be due to the reverse shock emission, 
if inverse Compton up-scatters primarily synchrotron optical photons into the X-ray 
range. The decay starts after the reverse shock has crossed the ejecta and electrons 
are no longer accelerated, and may have both a line of sight and an off-axis component 



CONTENTS 



41 



234 . This poses strong constraints on the Compton-y parameter, and cannot explain 
decays much steeper than a = —2, or —2 — /3 if the off-axis contribution dominates. 
Models involving bullets, whose origin, acceleration and survivability is unexplained, 
could give a prompt decay index a ~ — 3 to —5 [7"""|, with a bremsstrahlung energy 
index (3 ~ which is not observed in the fast decay; switching to a synchrotron or 
IC mechanisms requires additional parameters. Finally, a patchy shell model, where 
the Lorentz factor is highly variable in angle, would produce emission with a ~ —2.5. 
Thus, such mechanisms may explain the more gradual decays, but not the more extreme 
a = —5, —7 values encountered in some cases. It should be noted, however, that the 
Swift X-ray observations suggest that the steep decay is a direct continuation of the 
prompt emission |34Uj . which in turn suggests that the prompt and the fast decaying 
emission arise from the same physical region, posing a problem for the models in this 
paragraph (but not for the high latitude emission model). 

6.2. Shallow decay 

The slow decay portion of the X-ray light curves (a ~ —0.3 — 0.7), ubiquitously detected 
by Swift, is not entirely new, having been detected in a few cases by BeppoSAX. 
This, as well as the appearance of wiggles and flares in the X-ray light curves after 
several hours were the motivation for the "refreshed shock" scenario |405[ 1437 ] ( ffi).3j) . 
Refreshed shocks can flatten the afterglow light curve for hours or days, even if the 
ejecta is all emitted promptly at t = T < t 7 , but with a range of Lorentz factors, say 
M(T) oc T _s , where the lower T shells arrive much later to the foremost fast shells which 
have already been decelerated. Thus, for an external medium of density p oc r~ k and a 
prompt injection where the Lorentz factor spread relative to ejecta mass and energy is 
M(T) oc r~ s , -E'(r) oc r~ s+1 , the forward shock flux temporal decay is given by 437 

a = \(k - 4)(1 + s) + /3(24 - 7k + s/e)]/[2(7 + a - 2k)} , (38) 

(for more details, see Table"]"]). It needs to be emphasized that in this model all the 
ejection can be prompt (e.g. over the duration ~ T of the gamma ray emission) but the 
low T portions arrive at (and refresh) the forward shock at late times, which can range 
from hours to days. I.e., it is not the central engine which is active late, but its effects 
are seen late. Fits of such refreshed shocks to observed shallow decay phases in Swift 
bursts [180 lead to a T distribution which is a broken power law, extending above and 
below a peak around ~ 45. 

There is a different version of refreshed shocks, which does envisage central engine 
activity extending for long periods of time, e.g. ^ day (in contrast to the < minutes 
engine activity in the model above). Such long-lived activity may be due to continued 
fall-back into the central black hole |513j or a magnetar wind |522t 18*3*1 I338j . One 
characteristic of both types of refreshed models is that after the refreshed shocks stop 
and the usual decay resumes, the flux level shows a step-up relative to the previous 
level, since new energy has been injected. 
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From current analyses, the refreshed shock model is generally able to explain the 
flatter temporal X-ray slopes seen by Swift, both when it is seen to join smoothly on the 
prompt emission (i.e. without an initial steep decay phase) or when seen after an initial 
steep decay. Questions remain concerning the interpretation of the fluence ratio in the 
shallow X-ray afterglow and the prompt gamma-ray emission, which can reach ^ 1 [340J. 
This requires a higher radiative efficiency in the prompt gamma-ray emission than in 
the X-ray afterglow. One might speculate that this might be achieved if the prompt 
outflow is Poynting-dominated, or if a more efficient afterglow hides more of its energy 
in other bands, e.g. in GeV, or IR. Alternatively |211[ 1182] . a previous mass ejection 
might have emptied a cavity into which the ejecta moves, leading to greater efficiency 
at later times (although this would not work above the cooling frequency, which from 
the spectrum appears to be required in about half the cases), or otherwise the energy 
fraction going into the electrons increases oc t 1 / 2 . Other possible ways of addressing 
this include the afterglow coming from off-axis directions |lU4j . and exploring plausible 
reasons for having underestimated in previous studies the energy of the ejecta [182J. 

6.3. X-ray flares 

Refreshed shocks can also explain some of the X-ray flares whose rise and decay slopes are 
not too steep. However, this model encounters difficulties with the very steep flares with 
rise or decay indices a ~ ±5, ±7, such as inferred from the giant flare of GRB 0500502b 
(HZ] around 300 s after the trigger. Also, the flux level increase in this flare is a factor 
~ 500 above the smooth afterglow before and after it, implying a comparable energy 
excess in the low versus high V material. An explanation based on inverse Compton 
scattering in the reverse shock |234j can explain a single flare at the beginning of the 
afterglow, with not too steep decay. For multiple flares, models invoking encountering 
a lumpy external medium have generic difficulties explaining steep rises and decays 
[323l l528| . although extremely dense, sharp-edged lumps, if they exist, might satisfy the 
steepness 99 . 

Currently the more widely considered model for the flares ascribes them to late 
central engine activity [528[ 13381 l360j. The strongest argument in favor of this is that 
the energy budget is more easily satisfied, and the fast rise/decay is straightforward to 
explain. In such a model the flare energy can be comparable to the prompt emission, the 
fast rise comes naturally from the short time variability leading to internal shocks (or 
to rapid reconnection), while the rapid decay may be due to the high latitude emission 
following the flare, with to reset to the beginning of each flare (see further discussion 
in |B26j). However, some flares are well modeled by refreshed forward shocks, while in 
others this is clearly ruled out and a central engine origin is better suited [516] . Aside 
from the phenomenological desirability based on energetics and timescales, a central 
engine origin is conceivable, within certain time ranges, based on numerical models of 
the burst origin in long bursts. These are interpreted as being due to core collapse of a 
massive stellar progenitor, where continued infall into fast rotating cores can continue for 
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a long time |513| . However, large flares with a fluence which is a sizable fraction of the 
prompt emission occurring hours later remain difficult to understand. It has been argued 
that gravitational instabilities in the infalling debris torus can lead to lumpy accretion 
|371j . Alternatively, if the accreting debris torus is dominated by MHD effects, magnetic 
instabilities can lead to extended, highly time variable accretion |384| . which may give 
rise to GRB X-ray flares ftEE] . 

6-4- Late steep decay and jet breaks 

The late steep decay decay phase (4) of §H1 is seen in a modest fraction of the Swift 
bursts, mainly in X-rays, and mainly but not exclusively in long bursts. The natural 
interpretation is that these are caused by the fact that the outflow is collimated into 
a jet break: when the decrease of the ejecta Lorentz factor leads to the light-cone 
angle becoming larger than the jet angular extent, Tj(t) <; l/9j (e.g. S J5.5|) . the light 
curve steepens achromatically. For the Swift bursts, it is noteworthy that this final 
steepening has been seen in less than ~ 10% of the afterglows followed, and then with 
reasonable confidence mainly in X-rays. The corresponding optical light curve breaks 
have been few, and not well constrained. The UVOT finds afterglows in only ~ 30% 
of the bursts, and ground-based optical/IR telescopes have yielded few continued late 
time light curves monitored. This is unlike the case with the ~ 20 Beppo-SAX bursts, 
for which an achromatic break was reported in the optical |135j . while in rarer cases 
there was an X-ray or radio break reported, which in a few cases appeared to occur at 
a different time than the optical break (e.g. l 36J). 

The relative paucity of optical breaks in Swift afterglows may be an observational 
selection effect due to the larger median redshift, and hence fainter and redder 
optical afterglow at the same observer epoch. At higher redshifts the break occurs 
later in the observer frame, which compounds a possible reluctance to commit large 
telescope time on more frequently reported bursts (roughly 2/week from Swift versus 
an earlier 2/month with Beppo-SAX). One can speculate that the apparent scarcity of 
detected light curve breaks might indicate that at higher redshifts the opening angle 
is intrinsically larger. However, continued monitoring of the X-ray light curves with 
both Swift and Chandra is resulting in a growing number of bursts with high quality 
late X-ray light curves showing in some cases a clear break, and others the absence of a 
break up to weeks (also in short bursts, e.g. I184j ). This is an evolving topic, with 
some indications that light curve breaks may not (or not always) appear achromatic 
fTTTTl rTTTTj . 

6. 5. Prompt optical flashes and high redshift afterglows 

Optical/UV afterglows have been detected with the Swift UVOT telescope in roughly 
half the bursts for which an X-ray afterglow was seen. For a more detailed discussion 
of the UVOT afterglow observations see |411j . Of particular interest is the ongoing 
discussion on whether "dark GRB" are really optically deficient, or the result of 
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observational bias [40j. Another puzzle is the report of a bimodal intrinsic brightness 
distribution in the rest-frame R-band j269| I333j . This suggests possibly the existence of 
two different classes of long bursts, or at least two different types of environments. 

Compared to a few years ago, a much larger role is being played by ground-based 
robotic optical follow-ups, due to the increased rate of several arc-second X-ray alerts 
from XRT, and the larger number of robotic telescopes brought on-line in the last years. 
For the most part, these detections have yielded optical decays in the ^ few 100 s range, 
initial brightness my ~ 14 — 17 and temporal decay slopes a ~ 1.1 — 1.7 previously 
associated with the evolution of a forward shock |131|. 13*9*] . In a few prompt 
optical detection was achieved in the first 12-25 s |428[ I429l "480j. 
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Figure 8. The X-ray afterglow of the GRB 050094 at z = 6.29 [HH, showing for 
comparison the flux level of one of the most lumnious X-ray quasars at a comparable 
redshift, SDSS J1030+524 (multiplied by 100). The inset shows the GRB variability 
in the 10-70 ks timeframe. 

The most exciting prompt robotic IR detection (and optical non-detection) is that 
of GRB 050904 gHEEE]- This object, at the unprecedented high redshift of z = 6.29 
[220J , has an X-ray brightness exceeding for a day that of the brightest X-ray quasars (see 
Figure EJ) 491J . Its O/IR brightness in the first 500 s (observer time) was comparable 
to that of the extremely bright {my ~ 9) optical flash in GRB 990123, with a similarly 
steep time-decay slope a ~ 3 [Slj. Such prompt, bright and steeply decaying optical 
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emission is expected from the reverse shock as it crosses the ejecta, marking the start 
of the afterglow [HJ ED Ell] ■ 

However, aside from the two glaring examples of 990123 and 050904, in the last 
six years there have been less than a score of other prompt optical flashes, typically 
with more modest initial brightnesses m v ^ 13. There are a number of possible reasons 
for this paucity of optically bright flashes, if ascribed to reverse shock emission. One 
is the absence or weakness of a reverse shock, e.g. if the ejecta is highly magnetized 
305 . A moderately magnetized ejecta is in fact favored for some prompt flashes |521j . 
Alternatively, the deceleration might occur in the thick-shell regime (T 3> td ec . see 
eq. (|3f)jl . which can result in the reverse shock being relativistic, boosting the optical 
reverse shock spectrum into the UV |231j (in this case a detection by UVOT might 
be expected, unless the decay is faster than the typical 100-200 s for UVOT slewing 
and integration). Another possibility, for a high comoving luminosity, is copious pair 
formation in the ejecta, causing the reverse shock spectrum to peak in the IR |312j. 
Since both GRB 990123 and GRB 050904 had E iso ~ 10 54 erg, among the top few 
percent of all bursts, the latter is a distinct possibility, compatible with the fact that 
the prompt flash in GRB 050904 was bright in the IR I-band but not in the optical. 
On the other hand, the redshift z = 6.29 of this burst, and a Ly-a cutoff at ~ 800 
nm would also ensure this (and GRB 990123, at z = 1.6, was detected in the V-band). 
However, the observations of optical flashes in these two objects but not in lower Ei SO 
objects appears compatible with having a relativistic (thick shell) reverse shock with 
pair formation. Even in the absence of pairs, more accurate calculations of the reverse 
shock |326t 290J find the emission to be significantly weaker than was estimated earlier. 
Another possibility is that the cooling frequency in reverse shock is typically not much 
larger than the optical band frequency. In this case the optical emission from the reverse 
shock drops to zero very rapidly soon after the reverse shock has crossed the ejecta and 
the cooling frequency drops below the optical and there are no electrons left to radiate 
in the optical band [290 . 

7. Short GRB in the Swift Era 

7.1. Short GRB observations 

Swift, and in smaller numbers HETE-2, have provided the first bona fide short burst 
X-ray afterglows followed up starting ~ 100 s after the trigger, leading to localizations 
and redshifts. In the first of these, GRB 050509b j!55j the extrapolation of the prompt 
BAT emission into the X-ray range, and the XRT light curve from 100 s to about 1000 
s (after which only upper limits exist, even with Chandra, due to the faintness of the 
burst) can be fitted with a single power law of a ~ 1.2 (1.12 to 1.29 90% conf), or 
separately as a bat = 1-34 (0.27 to 2.87 90% conf) and a XRT =1-1 (0.57 to 2.36 90% 
conf). The X-ray coverage was sparse due to orbital constraints, the number of X- 
ray photons being small, and no optical transient was identified, probably due to the 
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faintness of the source. An optical host was however identified, an elliptical galaxy 
[53]. The next one, discovered by HETE-2, was GRB 050709 @EE|. Its host [M]] is an 
irregular galaxy at z — 0.16 (and the observations ruled out any supernova association). 
Even earlier, HETE-2 reported the short GRB 040924 [469 , with a soft gamma-ray 
prompt emission and a faint broken power law optical afterglow [202*] . A proposed host 
galaxy at z — 0.86 shows star formation, and evidence for an associated 1998bw-like SN 
contribution to the light curve |448j . which suggests this is perhaps the short end of the 
long burst or XRF distribution. The next Swift short burst, GRB 050724, was relatively 
bright, and besides X-rays, it also yielded both a decaying optical and a radio afterglow 
[3*"""j . This burst, together with a significant part of other short bursts, is associated 
with an elliptical host galaxy. It also had a low-luminosity soft gamma-ray extension 
of the short hard gamma-ray component (which would have been missed by BATSE), 
and it had an interesting X-ray afterglow extending beyond 10 5 s [2E] (Figure H"}. The 
soft gamma-ray extension, lasting up to 200 s, when extrapolated to the X-ray range 
overlaps well with the beginning of the XRT afterglow, which between 100 and 300 s 
has a ~ —2, followed by a much steeper drop a ~ — 5 — 7 out to ~ 600s, then a more 
moderate decay a ~ — 1. An unexpected feature is a strong flare peaking at 5 x 10 4 s, 
whose energy is 10% of the prompt emission, while its amplitude is a 10 times increase 
over the preceding slow decay. Among more recent Swift short bursts, such as GRB 
050813 407] I443j had an X-ray afterglow, a possible elliptical host, and was reported to 
be near a galaxy cluster at z — 1.7 — 1.9 [HH]. GRB 051210 |253j had an X-ray power 
law afterglow, with bumps or flares, and optical identifications still under consideration. 
GRB 051221a had X-ray and optical afterglows , and the host is a star forming galaxy 
at z = 0.55 [33Sj. 

7.2. Short GRB prompt and afterglow emission 

The main challenges for an understanding of the mechanism of short bursts are the 
relatively long, soft tail of the prompt emission, and the strength and late occurrence 
of the X-ray bumps or flares. A possible explanation for the extended long soft tails 
(~ 100s) may be that the compact binary progenitor is a black hole - neutron star 
system [26 , for which analytical and numerical arguments ([S], and references therein) 
suggest that the disruption and swallowing by the black hole may lead to a complex 
and more extended accretion rate than for double neutron stars (c.f. |315[ 1109"] ). The 
flares, for which the simplest interpretation might be to ascribe them to refreshed shocks 
(compatible with a short engine duration T ^ t 7 ~ 2 s and a distribution of Lorentz 
factors), requires the energy in the slow material to be at least ten times as energetic 
as the fast material responsible for the prompt emission, for the GRB 050724 flare at 
10 4 s. The rise and decay times are moderate enough for this interpretation within the 
errors. On the other hand, if the decay slope is -2.8, this is steeper than expected for 
refreshed shocks, but consistent with the high-latitude —2 — (3 model; a time origin to 
can be determined at the beginning of the flare, and late Chandra observations indicate 
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that the decay after the resumes where it had left off before the flare, which is more 
consistent with a late engine activity interpretation |27Uj . requiring a factor 10 less 
energy budget than the refreshed shock interpretation. Another interpretation for such 
flares might be an accretion-induced collapse of a white dwarf in a binary, leading to 
a flare when the fireball created by the collapse hits the companion [282] . which might 
explain moderate energy one-time flares of duration < 10 2 s. However, for repeated, 
energetic flares, as also in the long bursts, the total energetics are easier to satisfy if one 
postulates late central engine activity (lasting at least half a day), containing ~ 10% of 
the prompt fluence A possible way to produce this might be temporary choking 
up of an MHD outflow [384] (c.f. [478] ). which might also imply a linear polarization of 
the X-ray flare [lllj . Such MHD effects could plausibly also explain the initial ~ 100 
s soft tail. Another magnetic mechanism proposed for late X-ray flares in short bursts 
invokes a temporary post-merger massive neutron star [HEj- However, a justification for 
substantial <; 10 5 s features remains so far on rather tentative grounds. 
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Figure 9. The afterglow of GRB 050724 j^U, showing the Swift results on the prompt 
BAT emission extrapolated to the X-ray range and the subsequent XRT emission, as 
well as the late Chandra follow-up. 

The similarity of the X-ray afterglow light curve with those of long bursts is, in 
itself, an argument in favor of the prevalent view that the afterglows of both long and 
short bursts can be described by the same paradigm, independently of any difference 
in the progenitors. This impression is reinforced by the fact that the X-ray light curve 



CONTENTS 



48 



temporal slope is, on average, that expected from the usual forward shock afterglow 
model, and that in GRB 050724 the X-ray afterglow shows what appears like an initial 
steep decay, a normal decay and a significant bump or flare. The identification of jet 
breaks in short bursts is still preliminary, and the subject of debate. In two short bursts 
(so far) evidence evidence has been reported for a jet break [HEl I350| IBTJj . (However, 
in GRB 050724 a late Chandra observation indicates no X-ray break |184j ). Taking 
these breaks as jet breaks, the average isotropic energy of these SHBs is a factor ~ 100 
smaller, while the average jet opening angle (based on the two breaks) is a factor ~ 2 
larger than those of typical long GRBs |131| 1350] . Using standard afterglow theory, 
the bulk Lorentz factor decay can be expressed through T(td) = 6.5(n / E^qY^H^^ 8 , 
where = (t/day), n Q is the external density in units of cm -3 , and E$q is the isotropic 
equivalent energy in units of 10 50 ergs. If the jet break occurs at T(tj) = 9J 1 , for a 
single-sided jet the jet opening angle and the total jet energy Ej are 

9 j = 9°{n /E^Hf d , E 3 = 7T0j£ ~ 10 49 ny 4 (£ 50 t^) 3 / 4 erg . (39) 

For the afterglows of GRB 050709 and GRB 050724, the standard afterglow expressions 
for the flux level as a function of time before and after the break lead to fits 350 
which are not completely determined, allowing for GRB050709 either a very low or a 
moderately low external density, and for GRB050724 a moderately low to large external 
density. The main uncertainty is in the jet break time, which is poorly sampled, and 
so far mainly in X-rays. A better determined case of an X-ray light curve break is that 
of GRB 051221a, where combined Swift XRT and Chandra observations indicate a late 
break at tj ~ 5 days, leading to an estimated 9j ~ 15 degrees [59J. This is similar to 
jet angles calculated numerically for compact merger scenarios by |214t 15]. It is worth 
noting, however, that there are some indications that light curve breaks may not (or 
not always) be achromatic |113| 1361] . We note that chromatic breaks have been argued 
for in some long bursts, e.g. GRB 030329, suggesting different beam opening angles for 
the optical/X-ray and the radio components ([HZ]; see also |S70j), and independently of 
whether this is the explanation, a similar phenomenon may be present in short bursts. 

7. 3. Short burst hosts and progenitors 

The most dramatic impact of Swift concerning short GRB, after the discovery 
and characterization of the afterglows, has been in providing the first significant 
identifications of host galaxies, with the implications and constraints that this puts 
on the progenitor issue. Out of ten short bursts detected until the end 2005, four of the 
hosts (GRB 040924, 050509b, 050724 and 050813; \$M 123 ESI) are elliptical galaxies, 
one (GRB 050709, [20!) is a nearby irregular galaxy, and one (GRB 050906, [2E!) is 
a star-forming galaxy. The number of elliptical hosts is of significant interest for the 
most frequently discussed progenitor of short GRB, the merger of neutron star binaries 
343| I105[ EHl 1261] , which would be relatively more abundant in old stellar population 
galaxies such as ellipticals. The argument partly depends on the expected long binary 
merger times, which in early population synthesis and merger simulations [1H] was 
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taken to be in excess of 10 8 years. More recent populations synthesis calculations [3U] 
have reduced this to the point where compact mergers could be expected in substantial 
numbers also in young, e.g. star-forming galaxies, although statistically most mergers 
would be expected in old galaxies. The preponderance of claimed elliptical hosts, where 
star formation is absent, argues against alternative short burst origins such as short- 
lived outflows from massive stellar collapses |477j . The lack of any observed supernova 
emission weeks after the burst |197[ H3Uj also argues against a massive stellar collapse 
(where a Ib/c supernova could be expected), and also against a gravitational collapse 
of C/O white dwarfs to neutron stars leading to a supernova la jTHj- 

An alternative interpretation of short bursts is that they may be the initial brief, 
hard spike seen in giant flares of soft gamma repeaters, or SGR [2061 1349] . SGRs must 
be young objects, due to the fast field decay rate, and the total energy in giant SGR 
flares detected so far is at least two orders of magnitude too small to explain the short 
burst fluxes detected at z <; 0.2. The lack of recent star formation activity in the 
four mentioned elliptical hosts also indicates that at least some short bursts cannot be 
ascribed to SGRs. A statistical analysis indicates that the fraction of short bursts which 
could be due to SGRs is less than 15% (or less than 40% at 95% confidence level) |329j. 
It is interesting that a correlation analysis of short bursts with X-ray selected galaxy 
clusters |163j gives a better than 2a angular cross-correlation with clusters up to z = 0.1, 
which compared to model predictions would indicate that most short bursts originate 
within ~ 270 Mpc. Any connection between alternative candidates and a possible third 
category of bursts, intermediate between short and long 190J 132 [)\ I29fi[ I2()()j remains so 
far unexplored. 

7-4- Short burst redshifts and progenitor lifetimes 

With over a half dozen reasonably well studied short bursts (as of end of 2005), their 
distribution in redshift space and among host galaxy types, including both ellipticals 
and spiral/irregulars |383j . is similar to that of other old population objects, and thus 
is compatible with neutron star binaries or black hole- neutron star binaries |328j . This 
progenitor identification, however, cannot be considered secure, so far. Nonetheless, the 
most striking thing about short hard burst (or SHBs) hosts is that it includes a number of 
ellipticals, with low star formation rate (SFR), e.g. 050909b, 050724; and even for those 
SHBs with star forming hosts, e.g. 050709, 051221a, the SFR is lower than the median 
SFR for the long GRB hosts jSHl- This confirms that they are a distinct population, as 
indicated also by their intrinsic spectral-temporal properties versus those of long bursts 
I2HZ1 CHI - Using the BATSE flux distribution and the observed redshifts, the SHB 
local rates are inferred to be at least ~ 10 Gpc~ 3 yr -1 |328| I187j without beaming 
corrections, and larger including beaming. The progenitor lifetimes lead to interesting 
constraints, e.g. the simple time delay distribution P(t) oc t~ l expected from galactic 
double neutron star systems appears in conflict with the low average redshift of SHBs 
|1531 1328j . although it is not ruled out |187j . This has led to inferring a typical lifetime 



CONTENTS 



50 



for the progenitors of ~ 6 Gyr, and the suggestion that they might be neutron star-black 
hole binaries, rather than double neutron stars. However if the redshift z ~ 1.8 for GRB 
050813 is correct, the lifetime of the progenitor would be constrained to ^ 10 3 Gyr [39J. 
On the other hand, consideration of the star formation history of both early and late 
type galaxies suggests that at least half of the SHB progenitors have lifetimes in excess 
of ~ 10 Gyr |530j . Population synthesis models of double compact binaries [21] indicate 
two populations, with short (10~ 2 — 0.2 Myr and long (10 2 — 10 4 Myr) merger times, with 
NS-NS and BH-NS binaries distributed roughly 1:1 and 4:1 between these two merger 
time ranges, in apparent agreement with current SHB redshift and host distributions 
between ellipticals and SFR galaxies. The origin of a fraction of double neutron stars 
in globular clusters |183j would help to explain short bursts which are offset from their 
host galaxy. 

8. Long GRB Progenitors in light of Swift 

8.1. Long GRB hosts and progenitors 

Out of the ~ 90 long bursts (i 7 ^ 2 s) detected by Swift up to the end of 2005, in all cases 
where a host galaxy was identified this was of an early type, usually a blue star-forming 
galaxy p)2"l I444(. 145 6j . This was also the case for the thirty-some cases measured by 
Beppo-SAX (e.g. |471j ) in the previous seven years. More recent observational studies 
have indicated also that the long GRB host galaxy metallicity is generally lower than 
that of the average massive star forming galaxies [262*1 1263] I456J. This has implications 
for the expected redshift distribution of GRB [334j (c.f. |5()9j ). indicating that ~ 40% of 
long GRB may be at z ^ 4. Long GRB may, in principle, be detectable up to z ^ 25 — 30 

1252 EE1 EES!. 

The preponderance of short-lived massive star formation in such young galaxies, as 
well as the identification of a SN Ib/c light curve peaking weeks after the burst in few 
cases, has provided strong support for a massive stellar collapse origin of long GRBs, 
as argued by |514t I346| I513j . The relatively long duration of the gamma-ray emission 
stage in these bursts (2 s ^ t 1 < 10 3 s) is generally ascribed to a correspondingly long 
duration for the accretion of the debris [2831 1380j falling into the central black hole 
which must form as the core of a massive star collapses. (For initial stellar masses in 
excess of about 28 — 3OM , the core is expected to collapse to a BH, e.g. |142j . while 
for smaller initial masses 1OM < M„ < 28M Q the collapsing core mass is below the 
Chandrasekhar mass and is expected to lead to a neutron star). This accretion onto the 
black hole feeds a relativistic jet, which breaks through the infalling core and the stellar 
envelope along the direction of the rotation axis. 

A related massive core collapse mechanism has been considered by |476| I477j taking 
into account MHD effects in the disk and BH, in which the basic accretion time is 
short enough to be identified with short bursts, but magnetic tension can result in 
suspended accretion leading to long bursts. A mechanism based on the shorting of 
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a charge separation built up around newly formed black holes has been discussed by 
|4211 1422j . Other mechanism invoked include collapse of a neutron star to a strange 
star (e.g. [6~5~ ( I3~" H llOlj . The most widely adopted scenario, from this list, is the first 
one, in which the long GRB derive their energy from either the gravitational energy 
liberated by the torus of debris accreting onto the central BH formed by the massive 
core collapse, or by the extraction of the rotational energy of the BH, mediated by the 
presence of the debris torus, whose accretion lifetime in both cases is identified with the 
duration of the "prompt" gamma-ray emitting phase of the burst. 

8.2. Supernova connection 

In the year following the launch of Swift no supernovae were identified in association 
with GRB. In fact, there are some upper limits on possible supernovae, the most 
notable ones being on the short bursts GRB 050509b [Mj and GRB 050709 flgfl . 
However, from the previous eight year period there are two well documented cases of 
supernovae associated with long bursts, and several more weaker cases, were the evidence 
suggests a long GRB-SN connection. The first evidence for a long GRB-supernova 
association was discovered in GRB 980425, which appeared associated with SN 1998bw 
|150| 1241] . This was a peculiar, more energetic than usual Type Ib/c supernova, where 
the apparently associated GRBs properties seemed the same as usual, except for the 
redshift being extremely small (z ~ 0.0085). This implied the lowest ever long GRB 
isotropic equivalent energy ~ 10 48 erg, which resulted in the association being 
treated cautiously. However, using SN 1998bw as a template, other possible associations 
were soon claimed through detection of reddened bumps in the optical afterglow light- 
curves after a time delay compatible with a supernova brightness rise-time, e.g. in 
GRB 980326(11, GRB 970228(1011 USB, GRB 000911 GRB 991208|S3|, GRB 

990712|1SB, GRB 011121jSOj, and GRB 020405 [332] . 

The first unambiguous GRB-SN association was identified in GRB 030329, at a 
redshift z = 0.169, through both a supernova light-curve reddened bump and, more 
convincingly, by measuring in it a supernova spectrum of type Ib/c (i.e. the same type 
as in 1998bw) 455} 1196] . As a corollary, this observation rules out the "supra-nova" 
mo del 14*8*5] . in which a core collapse to neutron star and a supernova was assumed to 
occur months before a second collapse of the NS into a BH and a GRB; the delay 
between GRB 030329 and SN 2003dh is less than two days, and is compatible with both 
events being simultaneous[196j . For pre-Swift GRB-SN associations, see, e.g. |519t I518j . 

More recently, Swift observed with all three instruments, BAT, XRT and UVOT, an 
unusually long (~ 2000 s), soft burst, GRB 060218 [HB, which was found to be associated 
with SN2006aj, a very nearby (z = 0.033) type Ic supernova [^I^I^ISraiSSl liTgi. 
This supernova light curve peaked earlier than most known supernovae, and its time 
origin can also be constrained to be within less than a day from the GRB trigger. This 
is the first time that a connected GRB and supernova event has been observed starting 
in the first ~ 100 s in X-rays and UV/Optical light, and the results are of great interest. 
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The early X-ray light curve shows a slow rise and plateau followed by a drop after 
~ 10 3 s, with a power law spectrum and an increasing black-body like component which 
dominates at the end. The most interesting interpretation involves shock break-out of 
a semi-relativistic component in a WR progenitor wind [HI] (c.f. |114t 1461] ). After this 
a more conventional X-ray power law decay follows, and a UV component peak at a 
later time can be interpreted as due to the slower supernova envelope shock. Another 
GRB/SN detection based on Swift afterglow observations is that associated with GRB 
050525A E3]. 

8. 3. Jet dynamics, cocoons and progenitors 

For both long and short bursts, the most widely discussed central engine invokes a central 
black hole and a surrounding torus, either produced by a massive stellar core collapse 
(long bursts) or the merger of NS-NS or NS-BH binaries (short bursts). The latter 
mechanism is observationally on a less firm footing than the first, and in both collapse 
and merger cases the black hole could be preceded by a temporary massive, highly 
magnetized neutron star. There are two ultimate energy sources: the gravitational 
binding energy of the torus and the spin energy of the black hole. A possible third 
is the magnetic energy stored during the collapse, which derives its energy from the 
other two. Two main ways have been discussed for extracting the accretion energy and 
black hole spin energy, namely a neutrino-driven wind flOSl I419| U20, 283J 11451 1260J, 
and the Blandford-Znajek 46 mechanism. Both mechanisms lead to an optically thick 
e^ jet or fireball, but the second is dynamically Poynting-dominated, i.e. dominated 
by strong magnetic fields threading the black holejHOHl I278| I474| I266j . Needless to say, 
identification of the content of the fireball and the mechanism of GRB prompt emission 
would shed light on the mechanism that powers the central engine. Hence the excitement 
following claims of a very large gamma-ray polarization in GRB 021206 [Z2] suggesting 
a strongly magnetized central engine. This observation has been challenged |4231 151 lj . 
A strong gamma-ray polarization could in principle be expected from a pure Poynting- 
flux dominated jet [282], or in a baryonic hydrodynamic jet with a globally organized 
magnetic field configuration |499~1 I179| I177j . A strong but less extreme magnetization 
of the jet is inferred from a combined reverse-forward shock emission analysis of GRB 
990123 jnTTHTTU] . 

In all models, an e ± ,7 fireball is expected as a result of the dissipation associated 
with the transient core collapse or merger event. The initial chaotic motions and 
shears also are expected to lead to build up significant magnetic stresses [465J. A 
combination of the relativistic lepton (e 1 * 1 ) and MHD fields up to ~ 10 15 Gauss can 
provide the driving stresses leading to a highly relativistic expansion with Tj ^> 1. The 
fireball is very likely also to involve some fraction of baryons, and uncertainties in this 
"baryon pollution" |344j remain difficult to quantify until 3D MHD calculations capable 
of addressing baryon entrainment become available. If the progenitor is a massive star, 
the expectation is that the fireball will likely be substantially collimated, due to the 
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transverse containing pressure of the stellar envelope, which, if fast-rotating, provides a 
natural fireball escape route along the centrifugally lightened rotation axis. 

The development of a jet and its Lorentz factor in a collapsar has been discussed 
analytically in |311| 15041 I288| I25b'j . The essence of the dynamics of the jet in a burst 
from a massive star is that as long as the central BH accretes, it injects along the 
rotation axis a relativistic jet, whose dimensionless entropy must be comparable to 
or larger than the final bulk Lorentz factor of the jet once it has emerged from the 
star, r] = (L/Mc 2 ) i> I\, <; 100. Even though such a jet is highly relativistic as it 
is injected, the overburden of the stellar core and envelope slow the jet head down to 
a sub-relativistic speed of advance, which gradually increases as the jet moves down 
the density gradient of the star. The difference between the injection and advance 
speed causes gas and energy spill-over into a transrelativistic cocoon of waste heat 
|311| 1288} I25fij surrounding the jet, which may be detectable |392[ 1368] . By the time it 
reaches the boundary of the He core (Rue ~ 10 11 cm) the jet head has reached a speed 
v j ~ c. This takes, in the star's frame, ~ 10 s, hence the central engine must continue 
injecting energy and momentum into the jet for at least this long. A very sharp drop in 
density is predicted by stellar models at this radius, beyond which a tenuous hydrogen 
envelope extends as a power law. In going down this sharp gradient, the jet head Lorentz 
factor shoots up to a value comparable to its final value, Tj > 100 ( jSTJl l5*U4"] ).Once 
the jet head is relativistic, it becomes ballistic, and it is no longer affected by whether 
the central engine energy injection continues or not. A constraint on the mass of the 
envelope is that the mass overburden within the jet solid angle must be less than the 
jet total energy divided by Tjc 2 Q288J). If the star has lost is H envelope, this condition 
is guaranteed, e.g. as in Wolf-Rayet type stars, where a stellar wind phase leads to 
envelope loss previous to the core collapse phase. WR stars are, in fact, thought to 
be the progenitors of type Ib/c supernovae, which is the only type so far seen in a few 
cases associated with GRB. A modest envelope, however, should still be compatible 
with a high Lorentz factor, which could be tested through detection of weak H lines in a 
GRB associated supernova (and may also be tested through TeV neutrino observations, 

The 2D development of a relativistic jet making its way out through a star have 
been calculated numerically by, e.g. [5, 529], while magnetically dominated jets are 
discussed by |508t I102| 1385] . Jets in compact mergers have calculated numerically 
by |214[ Ej. The relativistic numerical calculations of GRB jets are, so far, mainly 
hydrodynamic, and involve approximations about the energy and momentum injection 
at the lower boundary, the numerical difficulties in covering the entire dynamical range 
being extreme. The results [529J show that a jet of Tj ~ 100 can escape a star of 
radius comparable to a WR (R* ~ 10 11 cm). The angular structure of the jet is, as 
expected, one where the Lorentz factor and energy per solid angle tapers off towards 
the edges, where instabilities cause mixing with and drag by the stellar envelope walls. 
An analytical argument |256j shows that this tapering off can result in an energy profile 
Ej{Q) oc 9~ 2 . Such a jet profile is a possible interpretation |414| I524j of the observational 
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correlation between the isotropic equivalent jet energy and the jet break time derived 
from a sample of burst afterglows |135t I352| . 

9. Very High Energy Photons and Non-Electromagnetic Emission 

The highly relativistic nature of the outflows is inferred from and constrained by the 
observations of GeV photons, which indicate the need for bulk Lorentz factors of V > 10 2 
|118[ 1191] I24j. Such Lorentz factors result in synchrotron spectra which in the observer 
frame extend beyond 100 MeV, and inverse Compton (IC) scattering of such synchrotron 
photons leads to the expectation of GeV and TeV spectral components |304j . While & 18 
GeV photons have been observed (e.g. |205j ). TeV photons are likely to be degraded to 
lower energies by 77 pair production, either in the source itself, or (unless the GRB is 
at very low redshifts) in the intervening intergalactic medium f73\. 192] . 

Besides emitting in the currently studied sub-GeV electromagnetic channels, GRB 
are likely to be even more luminous in other channels, such as neutrinos, gravitational 
waves and cosmic rays. For instance, nucleons entrained in the fireball will have 
<; 100 GeV bulk kinetic energies in the observer frame, which can lead to inelastic 
collisions resulting in pions, muons, neutrinos and electrons as well as their anti-particles. 
The main targets for the relativistic baryons are other particles in the relativistic 
outflow and particles in the external, slower moving environment. The expected flux 
and spectrum of 1-30 GeV neutrinos and 7-rays resulting from pion decay due to 
interactions within the expanding plasma depends, e.g., on the neutron/proton ratio 
and on fireball inhomogeneities, while that due to interactions with the surrounding 
medium depends on the external gas density and its distribution; and both depend 
on the Lorentz factor. Massive progenitors offer denser targets for nuclear collisions 
and a larger photon density for pj and 77 interactions, leading to modification of the 
photon spectra. On the other hand GRB from NS-NS mergers would be characterized 
by neutron-rich outflows, leading to stronger 5-10 GeV neutrinos and photons from np 
collisions [T7J E2 I416j . Photo-pion signatures of > 100 GeV photons and 10 14 — 10 18 
eV neutrinos may be expected to be relatively stronger in massive (high soft photon 
density) progenitors. Knowing what fraction of GRB, if any, arise from NS mergers 
is vital for facilitating interferometric gravitational wave detections, e.g. with LIGO. 
And, conversely, detection with LIGO would provide important clues as to whether 
short bursts are NS-NS (or NS-BH) mergers, or whether massive stellar collapses are 
asymmetric enough to produce substantial gravitational wave emission and serve as a 
test of the relationship between long GRB and supernovae. 

The Fermi mechanism in shocks developing in the GRB outflow can also accelerate 
protons to observer- frame energies up to ~ 10 20 eV j494[ 1492] . Internal shocks leading 
to the observed 7-rays have a high comoving photon density and lead to py photopion 
production and to > 100 TeV neutrinos [501]. In external shocks due to deceleration 
by the external medium, the reverse shock moving into the ejecta can produce optical 
photons ( §5.2j) which result in photopion production and <; 10 19 eV neutrinos 502J. 
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Neutrinos in the TeV to EeV range may be easier to detect than those at ~ 10 GeV 
energies, due to their higher interaction cross section, with instruments currently under 
construction. Such neutrinos would serve as diagnostics of the presence of relativistic 
shocks, and as probes of the acceleration mechanism and the magnetic field strength. 
The flux and spectrum of <; 10 19 eV neutrinos depends on the density of the surrounding 
gas, while the J> 10 14 eV neutrinos depend on the fireball Lorentz factor. Hence, the 
detection of very high energy neutrinos would provide crucial constraints on the fireball 
parameters and GRB environment. 

9.1. UHE photons from GRB 

Ultra-high energy emission, in the range of GeV and harder, is expected from electron 
inverse Compton in external shocks [304J as well as from internal shocks [365] in the 
prompt phase. The combination of prompt MeV radiation from internal shocks and a 
more prolonged GeV IC component for external shocks 303J is a likely explanation for 
the delayed GeV emission seen in some GRB [205J . (An alternative invoking photomeson 
processes from ejecta protons impacting a nearby binary stellar companion is |218j ). The 
GeV photon emission from the long-term IC component in external afterglow shocks 
has been considered by jHHl I523| 195} I488| 1489] . The IC GeV photon component is 
likely to be significantly more important |523j than a possible proton synchrotron or 
electron synchrotron component at these energies. Another possible contributor at 
these energies may be ir° decay from pj interactions between shock-accelerated protons 
and MeV or other photons in the GRAB shock region |531 146 7\ I137j . However, under 
the conservative assumption that the relativistic proton energy does not exceed the 
energy in relativistic electrons or in 7-rays, and that the proton spectral index is -2.2 
instead of -2, both the proton synchrotron and the components can be shown to be 
substantially less important at GeV- TeV than the IC component [523J . Another GeV 
photon component is expected from the fact that in a baryonic GRB outflow neutrons 
are likely to be present, and when these decouple from the protons, before any shocks 
occur, pn inelastic collisions will lead to pions, including 7r°, resulting in UHE photons 
which cascade down to the GeV range [HH UH I416j . The final GeV spectrum results 
from a complex cascade, but a rough estimate indicates that 1-10 GeV flux should be 
detectable [T7] with GLAST [Ml for bursts at z < 0.1. 

In these models, due to the high photon densities implied by GRB models, 77 
absorption within the GRB emission region must be taken into account [221 1272[ [398 , 
364, I365j . One interesting result is that the observation of photons up to a certain 
energy, say 10-20 GeV with EGRET, puts a lower limit on the bulk Lorentz factor 
of the outflow, from the fact that the compactness parameter (optical depth to 77) 
is directly proportional to the comoving photon density, and both this as well as the 
energy of the photons depend on the bulk Lorentz factor. This has been used by [272 
to estimate lower limits on T < 300 — 600 for a number of specific bursts observed with 
EGRET. On the other hand, for GRB with T > 850, TeV photons can escape the source 
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|33S] . 

Long GRB have recently been shown to be associated with supernovae ( A8.2J1 . If 
GRB also accelerate cosmic rays, as suspected, then these could leave long-lasting UHE 
photon signatures in supernova remnants which were associated with GRB at the time 
of their explosion. One example may be the SN remnant W49B, which may be a GRB 
remnant. A signature of a neutron admixture in the relativistic cosmic ray outflow would 
be a TeV gamma-ray signature due to inverse Compton interactions following neutron 
decay |2()9j (see also 13 ). Continued magnetic outflows upscattering companion 
photons may also signal GRB remnants |393j . The imaging of the surrounding emission 
could provide new constraints on the jet structure of the GRB. 

The recent detection of delayed X-ray flares during the afterglow phase of gamma- 
ray bursts (GRBs) with the Swift satellite (e.g. |528[ f338. 360J) suggests an inner-engine 
origin of these flares, at radii inside the deceleration radius characterizing the beginning 
of the forward shock afterglow emission. Given the observed temporal overlapping 
between the flares and afterglows, one expects an inverse Compton (IC) emission arising 
from such flare photons scattered by forward shock afterglow electrons |49Uj . The jet 
may also IC upscatter shock break-out X-ray photons 391J . This IC emission would 
produce GeV-TeV flares, which may be detected by GLAST and ground-based TeV 
telescopes. The detection of GeV-TeV flares combined with low energy observations 
may help to constrain the poorly known magnetic field in afterglow shocks. 

At higher energies, a tentative ^0.1 TeV detection at the 3a level of GRB970417a 
has been reported with the water Cherenkov detector Milagrito ^2]- Another possible 
TeV detection |37H| of GRB971110 has been reported with the GRAND array, at 
the 2.7a level. Stacking of data from the TIBET array for a large number of GRB 
time windows has led to an estimate of a ~ lo composite detection significance 
Better sensitivity is expected from the upgraded larger version of MILAGRO, as well as 
from atmospheric Cherenkov telescopes under construction such as VERITAS, HESS, 
MAGIC and CANGAROO-III IM21 HHJ EEHE1 B5BI 1123 1222- However, GRB 

detections in the TeV range are expected only for rare nearby events, since at this 
energy the mean free path against 77 absorption on the diffuse IR photon background 
is ~ few hundred Mpc [73| 192] . The mean free path is much larger at GeV energies, 
and based on the handful of GRB reported in this range with EGRET, several hundred 
should be detectable with large area space-based detectors such as GLAST [289, 523 . 

9.2. Cosmic rays from GRB 

In the standard fireball shock model of the prompt 7-ray emission, say from internal 
shocks or magnetic dissipation, and also in the external afterglow shocks, the same 
acceleration mechanisms which lead to the non-thermal electron power laws implied by 
the observed photon spectra must also lead to proton acceleration. Using the shock 
parameters inferred from broad-band photon spectral fits, one infers that protons can 
be accelerated to Lorentz factors up to ^ 10 11 in the observer frame |494[ 14821 . i.e. to 



CONTENTS 



57 



so-called GZK energy of E p ~ 10 20 eV. This is interesting mainly for "baryonic" jets, 
where the bulk of the energy is carried by baryons, whereas in Poynting-dominated jets 
there would be much fewer protons to accelerate. Well below the GZK energy, protons 
interacting with the MeV photons present in GRB or with thermal nucleons are above 
the pion production threshold and can produce ultra-high energy neutrinos, as discussed 
below. 

Discussions of GRB as cosmic ray sources are mainly oriented at exploring their 
contribution to the energy range above EeV (10 18 eV; e.g. [492 ), referred to as ultra- 
high energy cosmic rays, or UHECRs. (A model where GRB are responsible for CRs 
ranging from PeV to GZK is |512j ). At EeV and higher energies the observed UHECR 
isotropy and the small expected magnetic deflection suggests an extra-galactic origin. 
The requirement that they are not attenuated by the cosmic microwave background 
through photomeson interactions constrains that they are originated within a volume 
inside a radius of 50-100 Mpc, the so-called "GZK" volume (e.g. [IS]). Two broad 
classes of models suggested are the "top-down" scenarios, which attribute UHECR to 
decay of fossil Grand Unification defects, and the "bottom-up" scenarios, which assume 
UHECRs are accelerated in astrophysical sources. One of the most prominent candidate 
sources for bottom-up scenarios is GRBs [494l l482|l314| (two others are AGNs, e.g. [33] 
and cluster shocks, e.g. [208J). The most commonly discussed version of this scenario 
considers the UHECR to be protons accelerated in GRB internal shocks |494[ I493"] l492j. 
while another version attributes them to acceleration in external shocks |482t 1481} 199 j. 
(For UHECR acceleration in alternative GRB models, see, e.g. [8*% | 1 1 1 7 j ) . 

The persuasiveness of this scenario is largely based on two coincidences, namely, the 
required condition to accelerate protons to GZK energies is similar to the requirement 
for generating the prompt observed gamma-rays in GRB, and the observed UHECR 
energy injection rate into the universe (~ 3 x 10 44 erg Mpc -3 yr _1 ) is similar to 
the local GRB 7-ray energy injection rate [494, 1482j . These coincidences have been 
questioned, e.g. |459t I445j . but these objections have been resolved using new data and 
further considerations 492, 1481] . and GRBs remain a promising candidate for UHECRs. 
However, there are some caveats of principle. The internal shock scenario relies on the 
assumption that GRB prompt gamma-ray emission is due to internal shocks. Although 
this is the leading scenario, there is no strong proof so far, as is the case for the external 
shock (e.g., there are efficiency and spectrum issues, etc.). On the other hand, a Poynting 
flux dominated GRB model would have to rely on magnetic dissipation and reconnection, 
accelerating electrons and hence also accelerating protons- but details remain to be 
investigated. The external shock model would have to rely on a magnetized medium 
|481j to reach the desired cosmic ray energy (as expected in pulsar wind bubbles 236 
in the supranova scenario [485 , which however has become less likely since the almost 
simultaneous GRB 030329/SN 2003dh and the more recent GRB 060218/SN2006aj 
association). 

Direct confirmation of a GRB orgin of UHECRs will be difficult. The next 
generation cosmic ray detectors such as the Pierre Auger Observatory ^3] will have 
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a substantially enhanced effective target area, which will greatly improve the cosmic ray 
count statistics. This will help to disentangle the two scenarios (top-down or bottom-up) 
and will reveal whether a GZK feature indeed exists. Within the bottom-up scenario, 
the directional information may either prove or significantly constrain the alternative 
AGN scenario, and may eventually shed light on whether GRBs are indeed the sources 
of UHECRs. 

9.3. UHE neutrinos contemporary with gamma-rays 

Internal shocks in the GRB jet take place at a radius T{ ~ 2Tfc5t ~ 5 x 10 12 £t_3r3 00 
cm. Here Tj = 300 r 30 o is the bulk Lorentz factor of the GRB fireball ejecta and 
5t = 10~ 3 5t_3 s is the variability time scale. Observed 7-rays are emitted from the 
GRB fireball when it becomes optically thin at a radius <S ?Y Shock accelerated 
protons interact dominantly with observed synchrotron photons with ~MeV peak 
energy in the fireball to produce a Delta resonance, p'-f — > A + |501| . The threshold 
condition to produce a A + is E p E 7 = 0.2r 2 GeV 2 in the observer frame, which 
corresponds to a proton energy of E p = 1.8 x 10 7 i?~J kIeV r| 00 GeV. The subsequent 
decays A + — > rm + — > ra/i + z/ M — > ne + v e i)^v tl produce high energy neutrinos in the GRB 
fireball contemporaneous with 7-rays |501[ I388j . Assuming that the secondary pions 
receive 20% of the proton energy per interaction and each secondary lepton shares 1/4 
of the pion energy, each flavor of neutrino is emitted with 5% of the proton energy, 
dominantly in the PeV range. 

The diffuse muon neutrino flux from GRB internal shocks due to proton acceleration 
and subsequent photopion losses is shown as the short dashed line in Fig. E3 The flux 
is compared to the Waxman-Bahcall limit of cosmic neutrinos, which is derived from 
the observed cosmic ray flux [502 . The fluxes of all neutrino flavors are expected to be 
equal after oscillation in vacuum over astrophysical distances. 

The GRB afterglow arises as the jet fireball ejecta runs into the ambient inter- 
stellar medium (ISM), driving a blast wave ahead into it and a reverse shock back into 
the GRB jet ejecta. This (external) reverse shock takes place well beyond the internal 
shocks, at a radius r e ~ 4T 2 e cAt ~ 2 x 10 17 r| 50 At 3 o cm [HE]. Here T e « 250r 250 is the 
bulk Lorentz factor of the ejecta after the partial energy loss incurred in the internal 
shocks and At = 30At 30 s is the duration of the GRB jet. Neutrinos are produced in 
the external reverse shock due to pry interactions of internal shock accelerated protons 
predominantly with synchrotron soft x-ray photons produced by the reverse shock. The 
efficiency of pion conversion from p'j interactions in this afterglow scenario is much 
smaller than in the internal shocks |502j . 

In the case of a massive star progenitor the jet may be expanding into a wind, 
emitted by the progenitor prior to its collapse. In this case, the density of the 
surrounding medium, at the external shock radius, may be much higher than that 
typical ISM density of n ~ 1 cm -3 . For a wind with mass loss rate of 10 _5 M Q yr _1 
and velocity of v w = 10 3 km/s, the wind density at the typical external shock radius 
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Figure 10. Diffuse muon neutrino flux arriving simultaneously with the 7-rays 
from shocks outside the stellar surface in observed GRB (dark short-dashed curve), 
compared to the Waxman-Bahcall (WB) diffuse cosmic ray bound (light long-dashed 
curves) and the atmospheric neutrino flux (light short-dashed curves). Also shown 
is the diffuse muon neutrino precursor flux (solid lines) from sub-stellar jet shocks in 
two GRB progenitor models, with stellar radii ri.2.5 (H) and rn (He). These neutrinos 
arrive 10-100 s before the 7-rays from electromagnetically detected bursts (with similar 
curves for v^, v e and v T ) 396 . 

would be ~ 10 4 cm -3 . The higher density implies a lower Lorenz factor of the expanding 
plasma during the reverse shocks stage, and hence a larger fraction of proton energy 
lost to pion production. Protons of energy E p > 10 18 eV lose all their energy to pion 
production in this scenario 502J I484[ IRQ] producing EeV neutrinos. 

9. 4- Precursor neutrinos 

As discussed before, the core collapse of massive stars are the most likely candidates for 
long duration GRBs, which should lead to the formation of a relativistic jet initially 
buried inside the star. The jet burrows through the stellar material, and may or 
may not break through the stellar envelope [3 13^ Internal shocks in the jet, while it 
is burrowing through the stellar interior, may produce high energy neutrinos through 
proton-proton (pp) and photomeson (pry) interactions |396j . High energy neutrinos are 
produced through pion decays which are created both in pp and pj interactions. The 
jets which successfully penetrate through the stellar envelope result in GRBs (7-ray 
bright bursts) and the jets which choke inside the stars do not produce GRBs (7-ray 
dark bursts). However, in both cases high energy neutrinos produced in the internal 
shocks are emitted through the stellar envelope since they interact very weakly with 
matter. 

High energy neutrinos from the relativistic buried jets are emitted as precursors 
(~ 10-100 s prior) to the neutrinos emitted from the GRB fireball in case 
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of an electromagnetically observed burst. In the the case of a choked burst 
(electromagnetically undetectable) no direct detection of neutrinos from individual 
sources is possible. However the diffuse neutrino signal is boosted up in both scenarios. 
The diffuse neutrino flux from two progenitor star models are shown in Fig. EH one 
for a blue super-giant (labeled H) of radius R* = 3 x 10 12 cm and the other a Wof- 
Rayet type (labeled He) of radius R* = 10 11 cm. The Waxman-Bahcall diffuse cosmic 
ray bound [503], the atmospheric flux and the IceCube sensitivity to diffuse flux are 
also plotted for comparison. The neutrino component which is contemporaneous with 
the gamma-ray emission (i.e. which arrives after the precursor) is shown as the dark 
dashed curve, and is plotted assuming that protons lose all their energy to pions in py 
interactions in internal shocks. 

Most GRBs are located at cosmo logical distances (with redshift z ~ 1) and 
individual detection of them by km scale neutrino telescopes may not be possible. The 
diffuse v flux is then dominated by a few nearby bursts. The likeliest prospect for UHE 
v detection is from these nearby GRBs in correlation with electromagnetic detection. 
Detection of ultrahigh energy neutrinos which point back to their sources may establish 
GRBs as the sources of GZK cosmic rays. 

The detection of ultrahigh energy neutrinos by future experiments such as 
ICECUBE jUIl, ANITA [TT], KM3NeT |225| . and Auger [Hj can provide useful 
information, such as particle acceleration, radiation mechanism and magnetic field, 
about the sources and their progenitors. High energy neutrino astrophysics is an 
imminent prospect, with Amanda already providing useful limits on the diffuse flux 
from GRB HZ] and with ICECUBE HH EH2| on its way. The detection of TeV 
and higher energy neutrinos from GRB would be of great importance for understanding 
these sources, as well as the acceleration mechanisms involved. It could provide evidence 
for the hadronic vs. the MHD composition of the jets, and if observed, could serve as an 
unabsorbed probe of the highest redshift generation of star formation in the Universe. 

9.5. Gravitational waves 

The gamma-rays and the afterglows of GRB are thought to be produced at distances 
from the central engine where the plasma has become optically thin, r > 10 13 cm, 
which is much larger than the Schwarzschild radius of a stellar mass black hole (or 
of a neutron star). Hence we have only very indirect information about the inner 
parts of the central engine where the energy is generated. However, in any stellar 
progenitor model of GRB one expects that gravitational waves should be emitted 
from the immediate neighborhood of the central engine, and their observation should 
give valuable information about its identity. Therefore, it is of interest to study the 
gravitational wave emission from GRB associated with specific progenitors. Another 
reason for doing this is that the present and foreseeable sensitivity of gravitational 
wave detectors is such that for likely sources, including GRB, the detections would be 
difficult, and for this reason, much effort has been devoted to the development of data 
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analysis techniques that can reach deep into the detector noise. A coincidence between a 
gravitational wave signal and a gamma-ray signal would greatly enhance the statistical 
significance of the detection of the gravitational wave signal |125[ 1239 . It is therefore 
of interest to examine the gravitational wave signals expected from various specific 
GRB progenitors that have been recently discussed, and based on current astrophysical 
models, to consider the range of rates and strains expected in each case, for comparison 
with the LIGO sensitivity A general reference is [479J, which also discusses GRB- related 
sources of gravitational waves. 

Regardless of whether they are associated with GRBs, binary compact object 
mergers (NS-NS, NS-BH, BH-BH, BH-WD, BH-Helium star etc.) [Mil EZ21 EH E3H1 
E3 1213 Q3E1 HH and stellar core-collapses jSH3 El HI 12211 HSl WH\ have been 
studied as potential gravitational wave (GW) sources. These events are also leading 
candidates for being GRB progenitors, and a coincidence between a GW signal and a 
gamma-ray signal would greatly enhance the statistical significance of the former [T^o] . 
A binary coalescence process can be divided into three phases: in-spiral, merger, and 
ring-down |124[ I229j . For collapsars, a rapidly rotating core could lead to development 
of a bar and to fragmentation instabilities which would produce similar GW signals 
as in the binary merger scenarios, although a larger uncertainty is involved. The GW 
frequencies of various phases cover the 10 — 10 3 Hz band which is relevant for the Laser 
Interferometer Gravitational- wave Observatory (LIGO) |271j and other related detectors 
such as VIRGO gHZ], GEO600 PHEI and TAMA300 H2] . Because of the faint nature 
of the typical GW strain, only nearby sources (e.g. within ~ 200 Mpc for NS-NS 
and NS-BH mergers, and within ~ 30 Mpc for collapsars) [22l|] have strong enough 
signals to be detectable by LIGO-II. When event rates are taken into account [147[ 128]. 
order of magnitude estimates indicate that after one-year operation of the advanced 
LIGO, one event for the in-spiral chirp signal of the NS-NS or NS-BH merger, and 
probably one collapsar event (subject to uncertainties), would be detected[229 . Other 
binary merger scenarios such as BH-WD and BH-Helium star mergers are unlikely to 
be detectable |229j . and they are also unfavored as sources of GRBs according to other 
arguments [339J. 

A time-integrated GW luminosity of the order of a solar rest mass (^ 10 erg) is 
predicted from merging NS-NS and NS-BH models |239[ I42()[ 132 lj , while the luminosity 
from collapsar models is more model-dependent, but expected to be lower f |143t I318J; 
c.f. |47fij ). Specific estimates have been made of the GW strains from some of the 
most widely discussed current GRB progenitor stellar systems |229j . The expected 
detection rates of gravitational wave events with LIGO from compact binary mergers, 
in coincidence with GRBs, has been estimated by |125| I126j . If some fraction of GRBs 
are produced by DNS or NS-BH mergers, the gravitational wave chirp signal of the 
in-spiral phase should be detectable by the advanced LIGO within one year, associated 
with the GRB electromagnetic signal. One also expects signals from the black hole ring- 
down phase, as well as the possible contribution of a bar configuration from gravitational 
instability in the accretion disk following tidal disruption or infall in GRB scenarios. 
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The most promising GW-GRB candidates in terms of detections per year are the 
DNS and BH-NS mergers j229| (Fig. [TT|) . based on assumed mean distances from the 
formation rates estimated by |147j . More recent rate estimates are in |479j . and rates 
incorporating new information relating to Swift short GRB detections are in |2*9*1 H328J . 
Other binary progenitor scenarios, such as black hole - Helium star and black hole 
- white dwarf merger GRB progenitors are unlikely to be detectable, due to the low 
estimates obtained for the maximum non-axisymmetrical perturbations. 




f [Hz] f [Hz] 

Figure 11. Gravitational wave strain from a double neutron merger (left) and a 
collapsar (right) compared to advanced LIGO sensitivity |229| . 

For the massive rotating stellar collapse (collapsar) scenario of GRB, the non- 
axisymmetrical perturbations are very uncertain, but may be strong |90"t ll44[ I476j . and 
the estimated formation rates are much higher than for other progenitors |144| IHOj. with 
typically lower mean distances to the Earth. For such long GRB the rate estimates must 
incorporate the beaming correction |479j. This type of scenario is of special interest, 
since it has the most observational support from GRB afterglow observations. For 
collapsars, in the absence of detailed numerical 3D calculations specifically aimed at 
GRB progenitors, estimates were made j22P] of the strongest signals that might be 
expected in the case of bar instabilities occurring in the accretion disk around the 
resulting black hole, and in the maximal version of the recently proposed fragmentation 
scenario of the infalling cores. Although the waveforms of the gravitational waves 
produced in the break-up, merger and/or bar instability phase of collapsars are not 
known, a cross-correlation technique can be used making use of two co-aligned detectors. 
Under these assumptions, collapsar GRB models would be expected to be marginally 
detectable as gravitational wave sources by the advanced LIGO within one year of 
observations. Figure ^2 depicts the characteristic GW strains for the double neutrons 
star merger and the collapsar model. 

Other calculations of massive stellar collapse GRB [476, 477] take into account 
MHD effects in the disk and BH. More general studies of massive stellar core collapse 
event gravitational wave emission are presented in |146j . considering both core collapse 
SN and the progenitors of long GRB. 

In the case of binaries the matched filtering technique can be used, while for sources 
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such as collapsars, where the wave forms are uncertain, the simultaneous detection 
by two elements of a gravitational wave interferometer, coupled with electromagnetic 
simultaneous detection, provides a possible detection technique. Specific detection 
estimates have been made |229| I479j for both the compact binary scenarios and the 
collapsar scenarios, 

Both the compact merger and the collapsar models have in common a high angular 
rotation rate, and observations provide evidence for jet collimation of the photon 
emission, with properties depending on the polar angle, which may also be of relevance 
for X-ray flashes. Calculations have been made |23Uj of the gravitational wave emission 
and its polarization as a function of angle expected from such sources. The GRB 
progenitors emit I — m — 2 gravitational waves, which are circularly polarized on the 
polar axis, while the + polarization dominates on the equatorial plane. Recent GRB 
studies suggest that the wide variation in the apparent luminosity of GRBs are caused 
by differences in the viewing angle, or possibly also in the jet opening angle. Since 
GRB jets are launched along the polar axis of GRB progenitors, correlations among the 
apparent luminosity of GRBs (L^(9) oc 9~ 2 and the amplitude as well as the degree of 
linear polarization P degree of the gravitational waves are expected, P oc 9 4 oc L~ 2 . At 
a viewing angle larger than the jet opening angle 9j the GRB 7-ray emission may not be 
detected. However, in such cases an "orphan" (see, e.g. |203t I53l(l395| ) long- wavelength 
afterglow could be observed, which would be preceded by a pulse of gravitational waves 
with a significant linearly polarized component. As the jet slows down and reaches 
7 ~ 9J 1 , the jet begins to expand laterally, and its electromagnetic radiation begins to 
be observable over increasingly wider viewing angles. Since the opening angle increases 
as ~ 7 _1 oc t 1//2 , at a viewing angle 9 > 9j, the orphan afterglow begins to be observed 
(or peaks) at a time t p oc 9 2 after the detection of the gravitational wave burst. The 
polarization degree and the peak time should be correlated as P oc t 2 . 

Gravitational wave burst searches are underway with LIGO. The results from the 
third science run pP searched for sub-second bursts in the frequency range 100-1100 Hz 
for which no waveform model is assumed, with a sensitivity in terms of the root-sum- 
square (rss) strain amplitude of h rss ~ 10~ 20 Hz _1//2 . No gravitational-wave signals were 
detected in the eight days of analyzed data for this run. The search continues, as LIGO 
continues to be upgraded towards it ultimate target sensitivity. 
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